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ABSTRACT 


In recent times, air pollution has considerably received increased attention all over 
the world. Sulphur dioxide (SO2) is one of the major atmospheric pollutants resulting from the 
combustion of fossil fuels. The adverse effects of acid rain due to emission of SO2 on human, 
vegetation and buildings are well known. With stringent environmental regulations, the 
removal of SO2 by dry scrubbing methods has come into prominence because of its superior 
performance over the traditional wet scrubbing methods, especially in concentration levels of 
parts per million (ppm) and sub-ppm. 

In the present work, a combined experimental and theoretical study has been carried 
out on the removal of SO2 using solid sorbents. Both types of sorbents have been tested for 
their sorption behavior: non regenerative such as Ca(OH)2, CaCOa and Ca(CH3COOH)2, and 
regenerative such as zeolites. Reactions have been carried out in a packed bed of sorbent 
materials under varying operating conditions such as temperature, flowrate, and 
concentration, etc. 

CaO obtained on calcination of calcium acetate was found to be superior amongst all 
sorbents/adsorbents used in this work, with relatively longer breakthrough time and higher 
sulphate conversion ( 70 %). Although 5 A zeolites exhibited much shorter breakthrough time 
relative to calcium based sorbents under identical operating conditions, the former may still 
be preferred to non-regenerative types of sorbents due to their inherent regenerative 
characteristics. Use of a two-bed system, similar to PSA (pressure swing adsorption), may 
also be advantageous in the case of zeolite adsorbents. 

Mathematical model has been developed to understand the mechanism of sorption by 
sorbent materials having bi-disperse pores (as in zeolites). The theoretical analysis 
incorporates the effects of gas diffusion in the adsorbent bed, in the macro-pore and within the 
crystal. A simplified approach based on linear driving approximation and average radial 
concentrations within the macro and micro-pores is used. The model predictions are observed 
to be in good agreement with the experimental data. 
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CHAPTER 1 
INTRODUCTION 


Atmospheric emission of sulphur dioxide (SO 2 ) is primarily due to the combustion 
of fossil fuels such as coal and oil, which contain appreciable quantities of sulphur. Electric 
power plants are the major industrial users of these fossil ftiels [Liu and Liptak, 1980], On 
a worldwide perspective, SO 2 emissions are the dominant precursors of acid rain formation. 
Figurel.l schematically summarizes typical fate of a SO 2 species in atmosphere, beginning 
with the emission (anthropogenic and natural sources), transport to atmosphere and 
chemical and physical transformation leading to acid rain and the formation of aerosols in 
atmosphere. 

There are at present two major widely accepted technologies to control SO 2 
emission: wet scrubbing using suitable solvent and dry scrubbing using solid sorbents. Wet 
methods are generally advantageous over dry methods in terms of capital and operation 
cost and for the requirement of secondary treatment of wastes produced during the wet 
method. These processes however suffer from the disadvantages that they are not efficient 
if SO 2 concentrations are in low levels i.e. parts per million (ppm) and sub ppm. If 
operating parameters are optimized to improve the removal efficiency for low levels of SO 2 
concentrations in the vent gas, the process often becomes non-economical. Further, these 
processes are generally not energy efficient. As a matter of fact many industries employing 
wet scrubbing techniques had to be shut down as they failed to keep the emissions within 
environmental regulations [Flood, 1971]. 

Dry scrubbing methods using suitable sorbents have been in practice for some time 
for controlling emission of acid gases such as SO 2 . Most of the earlier works done in this 
area, however, focussed on the aspect of corrosion on the interior sidewalls of the boiler, 
which recovered waste heat from gases emitted from the coal-based power plants. Among 
these works Boll and Patel (1961) analyzed the problem of corrosion and gas-side fouling in 
turbines and boilers by determining the composition of the flue gases. On the same line, 
Cutler et al. (1980) have analyzed the problem of corrosive gases environment created in 
fire-side boiler caused by combustion gases and ash particles. With increasing stringent 



environmental regulations some studied were carried out only in recent years to develop 
suitable technology to control SO 2 emissions in low levels. 

According to the present environmental regulations the atmospheric concentration 
of SO 2 emitted by the industries is not allowed to exceed 0.3 ppm on an annual mean 
concentration basis [Seinfeld and Pandis, 1998]. Dry scrubbing usually employs a packed 
bed reactor of solid sorbents for the reduction of these emissions. Both non-regenerative as 
well as regenerative types of materials may be used. Common examples of regenerative 
type of sorbents are calcium-based sorbents such as calcium carbonate, calcium hydroxide, 
and calcium acetate. The sorption process is usually carried out at high temperatures. The 
sorbents, however, suffer from rapid loss of reactivity and incomplete utilization. In the 
present study sulfation of these commercially available materials following calcination at 
high temperatures has been carried out. Precipitated CaCOs is prepared by calcining CaC 03 
followed by making Calcium hydroxide solution in water and passing CO 2 gas in the 
solution. The effect of CaCl 2 as an additive in precipitate CaO on SO 2 capturing efficiency 
is also studied. The effects of various operating variables such as types of sorbents, reaction 
temperature, SO 2 inlet gas concentrations, sorbent particle size on SO 2 capturing efficiency 
are ascertained. 

In recent times, zeolites have drawn considerable interest in air pollution 
applications. The studies have shown that a number of air pollutants such as SO 2 , NOx and 
volatile organic compounds (VOCs) may be removed using zeolites [Kopac et al. (1997), 
Srinivasan et al. (1999), Coker et al. (1999)]. The major advantage of using such type of 
materials is their ability for regeneration. The reaction temperature is also low. In the 
present study, SO 2 sorption by 5A commercial zeolites has been carried out. The 
experimental variables are reaction temperature, SO 2 inlet gas concentrations and particle 
size. Mathematical models are developed primarily to predict breakthrough characteristics 
of the sorbent materials. The theoretical analysis takes into account intra-particle diffusion 
within the mono-pores of CaO based sorbents and bi-pores of zeolites materials. The 
removal process is mainly due to irreversible reaction in the case of CaO sorbents and 
reversible adsorption in the case of zeolites adsorbents. 


1.1 Objectives: 

The main objectives of the present study are as follows; ' 

1 ) design, fabrication and set-up of an experimental test bench to carry out the 
experimental study. 

2) determining optimal calcination temperature for Ca-based sorbents. 

3) determining breakthrough characteristic of the materials under varying operating 
conditions such as particle diameter, flue gas concentration, amount of sorbent and flow 
rate of the gas, on the removal of SO 2 . 

4) preparing the modified CaCOa by commercial CaCOj and study the breakthrough 
characteristic as same defined in steps 2 and 3. 

5) study the effect of CaCh on modified CaC 03 on the sorption of SO 2 removal process 

6) determining the optimal adsorption temperature for Zeolite . 

7) determinig the breakthrough characteristic of Zeolite on the adsorption of SO 2 by study 
the same operating condition as defined in steps 2 and 3 . 

8) development of a theoretical analysis and mathematical model for both types of 
adsorbent. 

9) model verification with experimental results for both types of sorbents 

1 0) a detailed breakthrough as well as oxide to sulphate conversion analysis of all the 
experiments . 

1 1 ) surface analysis (SEM and BET area measurements). 




Figurel.l: SO2: Sources , effects and control techniques 
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CHAPTER 2 

LITERATURE SURVEY 

In this chapter, a review of the experimental studies on SO 2 removal by non-regenerative 
(i.e. CaO) and regenerative (i.e. zeolites and activated carbon) types of sorbents has been 
presented. A review of various theoretical analysis carried out to understand the mechanism 
of sorption (in case of CaO) and adsorption (in case of zeolites) of SO 2 has also been 
presented. The studies done so far range from experimental investigations on change in 
surface morphology (pores structure) of the sorbents during reaction to ascertaining relative 
reactivity of various types of sorbents with SO 2 to preparation of sorbents of relatively larger 
surface area. Various mechanisms have also been proposed in these studies to explain the 
gas-solid reactions. An attempt has been made in the present work to identify areas, which 
require further study. 

Dogu et al. (1980) have studied the pore structure variation during the gas-solid 
reaction. He reported that the initial pore stmcture of the solid reactant has significant effects 
on the kinetics of reactions. The initial pore size distribution of solid reactant (CaO) was 
found to influence both the diffusion resistance of the gaseous reactant (SO 2 ) through the 
pores and the active surface area of the solid. The pore size distribution of the calcined CaO 
sample was also found to be affected considerably by the calcination temperature. The 
authors experimentally determined that the sample calcined at 950 "C had the lowest 
diffusion resistance and the highest activity for SO 2 removal. 

Dogu et al. (1982) studied the breakthrough of SO 2 during reaction with calcined 
limestone. In the experiment, a sigmoidal type of output SO 2 concentration vs. time profiles 
was observed. It was also observed that reaction temperature shifted the sigmoidal 
breakthrough curves to the left. The longest breakthrough time reported was at 720 °C and 
the corresponding maximum sulfate conversion was 35 %. 

Jorgensen et al. (1987) studied the reaction between sulfur dioxide and alkaline 
sorbents carried out in a fixed bed reactor at low temperature. Three m^or groups of 
sorbents were tested: clacitic and dolomitic hydrated lime, sodium bicarbonate, and lime 
stone (calcite and aragonite) mixed with salt additives such as calcium chloride (CaCb), 
sodium chloride (NaCl), and nitrates of calcium and sodium. Three major variables were | 
found to influence the SO 2 removal: SO 2 concentration in the flue gas, the dry bulb 
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temperature of the gas and relative humidity. The most notable result was a significant 
increase in SO 2 removal using CaCl 2 as an additive in an optimum proportion (O.lg/lg). It 
was believed that the CaCl 2 having a high water affinity (hygroscopic) increased the amount 
of water absorbed by the particle, thereby reducing the SO 2 ash blinding effect. Relative 
humidity was found to be critical in the removal of SO 2 . Lowering of the relative humidity 
from 70%to 50% stopped virtually all reactions of S 02 with CaCOa. 

Bruce et al. (1 989) experimentally compared the SO 2 reactivity of CaO obtained 
from CaCOa as well as Ca(OH) 2 . The variations of porosity and conversion during the 
course of the reaction were determined and conclusions were drawn that CaO obtained from 
Ca(OH )2 gave a superior performance than that from limestone. 

Yiannis et al. (1993) have studied the effectiveness of commercially available 
calcium magnesium acetate from cryogenic technology (Aldrich) (CMA) as sulphur 
capturing agent and explored the possibility of using CMA as a “ dry scrubbing” medium in 
boiler injection. In addition to versatility of the water-soluble CMA for being relatively 
easily sprayed during dry injection and therefore, eliminating grinding costs, there was 
additional useful property associated with CMA. CMA exhibited ability to form highly 
cenospheric, “popcom”-like oxide particles on heating at high temperature (750T150“C). 
These cenospheres posses thin porous walls with blowholes that enable penetration of SO 2 
in the interior of the particle and promotes complete utilization of the sorbent particles. The 
utilization of CM A depended on the cenospher wall thickness rather than the particle size 
and thus, outperformed other sorbents regardless of the size of resulting oxide particles. 

Mahuli et al. (1997) developed a new method to prepare CaCOa, which on calcination 
produced CaO of relatively larger surface area resulting in the surfate conversion of 70-75% 
compared to 30-35% for commercial sorbents under identical conditions. The sulfation 
characteristics were investigated at high temperatures: 900, 1000, and 1080°C and short 
contact times (20-600 ms) using small particle sizes (< 5|im) in entrained flow reactor. The 
modified carbonate was prepared by carbonation-precipitation in a slurry bubble column 
reactor using different surfactants in small concentration of about 2% by weight. The 
carbonate used in this study was prepared by bubbling CO 2 at a flowrate of 2.25 slpm for 20 
min through a Ca(OH )2 suspension (2.25wt.% solid concentration). The post-reaction sluny 
was filtered and dried in a vacuum oven at 75°C for 24 hr. The calcination and sulfation 
studies were performed at a fixed reaction temperature. Calcination was carried out under 
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inert N 2 environment. The gas composition for the sulfation runs was maintained at 5.45% 
O 2 , 3900 ppmv SO 2 and balance N 2 . The sulfation characteristic of this modified carbonate 
(MC) was compared with the commercial Linwood carbonate (LC) and a modified 
hydroxide (MH). It was observed that in spite of the pore properties of MH and MC being 
similar, the modified carbonate reactivity was much higher than that of MH. This was 
attributed to the evolution of pore structure. It was observed that, under conditions of 
interest an optimum pore-size range existed between 50-200 which provided sufficient 
surface area for the sulfation reaction without causing rapid pore filling and pore-mouth 
plugging. The development of ppt-CaO having larger surface area has also been reported 
elsewhere [Ghosh-Ddostidar et al., 1996]. 

In a recent work, Dasgupta et al. (2001) has carried out a combined experimental and 
theoretical study to understand the kinetics and mechanism of SO 2 removal by CaO 
sorbents. The breakthrough characteristics and sulphate conversion analysis were carried 
out for the sorption process. Mathematical models were developed for non-porous as well as 
porous materials. It was found that the model assuming a porous structure of CaO explained 
the breakthrough characteristics, even if the sorbent materials had small surface area (~1-10 
m^/g). It was observed that sorption process was best suited for low concentration levels of 
SO 2 i.e. in ppm and sub ppm. The favorable operating temperature for S02-CaO reactions 
was determined to be 950"C. Shorter breakthrough time and lower sulphate conversion at 
above 1000‘’C were attributed to breaking of the pore structure of the sorbent. No 
appreciable reaction was observed below 750“C. It was also observed that the particle of 
smaller diameter resulted in longer breakthrough and higher conversion. However, the 
pressure drop and particle agglomeration in the packed bed increased with decrease in the 
particle size. 

A literature survey shows that not much work has been done on the application of 
zeolites in control of SO 2 . Only in the past few years, some of the works have, beginning to 
emerge in this area, as environmental laws become stringent. 

Kopac et al. (1997) have studied adsorption behavior of SO 2 on commercially 
available molecular sieves (5A, 4A and AW300 Zeolites). The experiments were carried out 
in the temperature range of 523-718 K. The adsorption equilibrium parameters were 
obtained from non-isobaric pulse chromatographic experiments and using a method of 
“moment-analysis”. It was found that the equilibrium or partition constant decreased 
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considerably with increase in the temperature. Fro example, in the case of SO 2 adsorption on 
5 A molecular sieve, an increase of temperature from 523K to 673K resulted in decrease of 
the equilibrium constant from approximately 150 to 12. It was also found that the adsorption 
of SO 2 by various adsorbents used were fbund to decrease in the order of AW300 > 4A >5 A. 
The heats of adsorption of SO 2 on these zeolites determined experimentally were found to 
be in the order of 1 0.0 kcal/mol. 

In another study, Kopac et al. (1998) studied the adsorption behavior of SO 2 on 13X 
molecular sieve and activated carbon in a temperature range of 353-453K and 523-673K, 
respectively. The exothermic heats of adsorption of SO 2 on activated carbon and 13X 
molecular sieve 13X were determined as 8.99 and 12.4 kcal mol"', respectively. 

Tsibranska et al. (2000) developed a mathematical model for adsorption of SO 2 on 
natural clinoptilolite. These materials are bi-porous in nature having macro and micro pores, 
similar to zeolites molecular sieves. The model accounted for macro-pore diffusion, micro- 
particles size distribution, concentration dependent micro-pores diffusivity and nonlinear 
adsorption equilibrium on the micro-particle surface. The external mass transfer in the film 
around the particles and axial dispersion in the bed were also taken into account. The kinetic 
and equilibrium parameters of the model were extracted from the experimental data on 
breakthrough. The nature of micro particle size distribution was independently 
experimentally determined from SEM analysis. A good agreement between the experimental 
data and model results was obtained for larger pellets where macro-pore diffusion control 
was of greater importance. On the other hand, in micro-particles where micro-pore diffusion 
is rate limiting, the model prediction was observed to be less satisfactory. 

In this work, experimental studies on non-regenerative sorbents carried out by 
Dasgupta (2001) have been extended to include other types of Ca-based sorbents. The 
especially prepared precipitated calcium carbonate sorbents have also been tested for their 
sorption performance. The focus of the present work is on the application of regenerative 
type of sorbents like zeolites in control of SO 2 , which have shown potential in adsorbing 
several air pollutant gases, such as SO 2 , NOx, and volatile organic compounds (VOCs). 
Extensive theoretical and experimental studies have been carried out in this work to 
determine the breakthrough characteristic of 5 A zeolite molecular-sieve materials during the 
adsorption of SO 2 . 
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CHAPTER 3 

THEORETICAL ANALYSIS 

3.1 Non-regenerative sorbent (CaO): 

In non-regenerative types of the sorbents (CaO), the removal of SO 2 takes place 
through two reaction steps, a) calcination of the sorbent materials (CaCOa) to produce 
CaO of higher surface area and porosity, b) reaction with SO 2 (sulfation) in the presence 
of O 2 to form the higher molar volume solid product, CaS 04 . These steps are represented 
by the following chemical reactions; 

a, Ca(OH )2 -> CaO + H 2 O 

CaCO, Ca0 + C02 

b. CaO + SO 2 +1/2 O 2 _> CaS 04 

'I'he commercially available limestone is generally non-porous and possesses a 
very low surface area. However, on calcinations it yields a porous CaO, which is 
responsible for SO 2 capture. CaO reacts with SO 2 in the presence of oxygen to form a 
solid product, CaS 04 . The overall rate of reaction of this process is mainly dependent on 
CaS 04 formation on CaO. The molar volume of CaS 04 is sufficiently large enough so as 
to usually plug the pores prior to complete utilization of the sorbent. Pore plugging and 
thus the loss of porosity of the particle are the dominant cause of premature termination 
and incomplete utilization of the sorbent. 

The mathematical models for SO 2 removal by the sorbent materials assuming non- 
porous and porous structure have been described in details by Dasgupta (2001). Figures 
3.1 and 3.2 are the schematic representations of the overall sorption process in a non- 
porous and porous particle, respectively. In the porous-model, a particle or grain of CaO is 
assumed to consist of spherical sub-grains of uniform size separated by pores through 
which the reacting gases diffuse. As the reaction proceeds, a shell of the reaction products 
is gradually formed on the surface of the sub-grains, which intum offers resistance to the 
diffusion of SO 2 within the pores. 

The details of the mathematical formulation can be obtained from the 
aforementioned work [Dasgupta, 2001]. In brief, the theoretical analysis for the model 
assuming non-porous structures consisted of a species balance in the packed bed of CaO, 
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in the product sold layers of CaS 04 and on the surface of the un-reacted CaO. The 
variation in the radius of the reaction front of CaO was also incorporated in the analysis. 
It was assumed that an isothermal plug flow existed in the bed. The pressure drop was 
assumed to be negligible and there was no variation in the fluid velocity along the length 
of the bed. 

The theoretical analysis assuming a porous structure incorporates a species balance 
within the pores or sub-grains of CaO (diffusion of gases within the pores), in addition to 
that in the packed bed and on the surface of the material. As pointed out earlier in the text 
(Chapter 2), the model assuming a porous structure could explain the breakthrough 
characteristics of CaO over the entire experimental conditions used in the work. 

3.2 Regenerative sorbent (zeolite): 

Zeolites are crystalline aluminosilicates of alkali or alkali earth elements such as 
sodium, potassium, and calcium, represented by the stoichiometry; 

Mx/„[(A102)x(Si02)y]zH20 

where, x is the integer with y/x equal to or greater than 1, n is the valence of cation M, and 
z is the number of water molecules in each unit cell called sodalite cage [Yang et all 997], 
More than 1 50 types of zeolite have been synthesized to-date. These are designated by a 
letter or groups of letters: type A, X, Y, ZSM, and so on, depending upon the stmcture of 
the crystals. When the units are linked through four member rings zeolite of type A is 
formed. In the present work, SO 2 breakthrough analysis has been carried out on 5 A 
zeolite. Type A Zeolite contains 24 tetrahedral, 12 AIO4 and 12 Si 04 . When fully 
hydrated, 27 water molecules are contained in the central cage or cavity of the unit cell, 
and in eight smaller sodalite cages. 

Zeolite solids are available in the pellets form, and these pellets are made by 
compressing zeolite crystals together, usually with a small percentage of binder to join the 
crystals together [Do et al.l998]. This results in a porous structures having macro and 
micro pores. Figure 3.3 is the schematic representation of such a biporous particle. The 
crystals are of the order of 0. 1 to 1 micron, and the zeolite pellets are of that of 1 mm. The 
void between the microparticle contributes to the micro-pores of the particle. These pores 
act as conduit to transport molecules from the surrounding in to the interior of the particle. 
Once inside the particle, molecules adsorb at the pore mouth of the micropores and thence 
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the adsorbed species diffuse into the interior of the crystal. Micro pores within the crystal 
provide the adsorption space to accommodate adsorbate molecules. The diffusion process 
in the macropore and the micropore follows the combination of the molecular and 
Knudsen mechanisms while the diffusion process inside the crystal follows an intra 
ciystalline diffusion mechanism. 

The mathematical model developed to explain the breakthrough of SO 2 in the 
packed bed consisting of zeolite pellets is based on the following assumptions: 

1 ) The temperature is uniform throughout the bed, pellet and crystals 

2) Negligible pressure drop along the length of the bed 
3 ) There is constant fluid velocity through out the bed 
4) Linear adsorption isotherm 

3.2, a. Mass Balance in the packed bed reactor: 

Since the SO 2 concentration level is in sub-ppm level, therefore SO 2 exothermic heat of 
absorption is negligible and isothermal condition can be assumed [Yang et al, 1997]. 
Consider an isothermal plug flow system in a packed bed of spherical pellets. By making a 
SO 2 mass balance in the packed bed of porosity Cb, the following equation is obtained; 


dt Sft az Gf, 

A B CD 


(3.2.1) 


Here, A = Transient term, B = Dispersion term, C = Convective term, D = Diffusion flux 
term at the outer surface of the pellet and represents diffusion of the gas from the bulk into 
the macropores of the zeolite pellets. 


In equation (3.2.1), 


fiz) = KJC,-C^ J 


(3.2.2) 


Cb = concentration of reactant gas at an arbitrary location z in the bed. 
Dz = axial dispersion co-efficient in the packed bed. 

V/.= superficial velocity 


a* = external surface area per unit volume of the pellet 


3(1-6,) 


Km = average mass transfer co-efficient 

Cp= concentration of the reactant gas in the macropores of the pellet 
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Rp = radius of the pellet 

z = axial distance along the bed 

The boundary conditions for the above equation are 

at t = 0, Cb = 0 

at t>0, z = 0 Cb=Cbo 


dz 


(3.2.3) 

(3.2.4) 

(3.2.5) 


3.2.b. Mass balance at any location R within the macro-pore volume of the pellet 
located at a distance of z in the bed: 

Considering radial diffusion within the pore volume and adsorption at the pore mouth of 
the crystals, the following equation is obtained: 


a 



E 


1 d(/i%) 
dR 
F 


+ 3(l-a)— = 0 
dt 


G 


(3.2.6) 


E = Transient term, F = Diffusion flux in the radial direction of the pellet, G = Rate of 
adsorption at the pore-mouth (or opening) of the crystals. 

Here, a = void fraction of the macro-pores 

Cp- gas concentration in the macropores of pellet 
Nr = molar diffusion flux in the radial direction of the pellet 
q concentration of the adsorbate inside the crystal 

q = average amount adsorbed at a radial location R within the pellet 

3.2.C. Mass balance at any location r within the micro-pore of the crystal located at a 
distance of R in the pellet: 

Mass balance in the crystal is given by the following solid diffiision equation; 
djr'^dqldr) 

dr C^^^) 


where Do= diffusivity of the adsorbate within the crystal 


The governing equations (3.2.1), (3.2.6) and (3.2.7) form the basis of the model developed 
in this study for predicting the adsorbents performance, especially breakthrough, in the 
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bed. As these equations are coupled partial differential equations (independent variables 
being time, axial and radial directions, both in r and R), an approach is adopted in the 
present work to simplify numerical computations and significantly reduce the CPU time. 
Essentially, in this approach radial (r and R) concentration profiles within the solid pores 
(both macro and micro) are averaged and then, the average gas phase concentrations 
within the pores of the pellets and the average adsorbate phase concentrations within the 
crystals are determined [Yang et al.l997, Gupta, 2002], The salient advantage of this 
mathematical approximation is the reduction of 2"*^ order PDEs to 1“^ order PDE with 
variation only in z direction. Moreover, the number of the equations is reduced to two. 
The detailed computational steps are as follows: 

First, the particle volume averaged quantities are calculated as follows: 



(3.2.8) 



where, the terms on the left hand side of the expressions are the volume average surface 
concentration in the pores and within the crystals (radius Rc) , respectively. 

Second, by integrating equation (3.2.6) with respect to R and using the volume averaged 
quantities the following equation is obtained: 

where, =Average concentration of reactant gas in the macropores of the pellet 


q = Averaged concentration of the adsorbate over the entire pellet and defined as 
volume average concentration: 

q = -^lqR^dR 

^ p 

The required initial and boundary conditions for Eq. (3.2.10) are as follows: 

At t = 0, ^ = 0 and^=0 


(3.2.12) 
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and t >0, 


at R = 0 


dR 


= 0 


(3.2.13) 


In the integration of Eq. (3.2.6), the following boundary conditions were used: = 0 at R 

= 0, due to symmetricity. At the surface (R = Rp), the diffusion flux into or out of the 

K-^R, is given by the boundary condition in pellet: 


pores, Np 
N 






(3.2.14) 


Substituting equation (3.2.14) in equation (3.2.10) we obtain 
(ic'’,, 3 ^ , 3(l-a) 9^ 

i T- - O 


(3.2.15) 


C = (’ + DR 


dt Rp a ' ■' a dt 

Third, parabolic concentration profiles for Cp and q were assumed within a pellet as 
follows: 

(3.2.16) 

q = A + (3.2.17) 

Fourth, this has been mathematically shown that a parabolic concentration profile q(r) 
combined with the solid diffusion equation (3.2.7) leads to the following linear driving 
force (LDF) approximated equation [Yang et al.l997]: 
dq _ 15/7. 


dt R 


-</) 


(3.2.18) 


Taking the volume-average of Eq. (3.2.18) over the entire pellet we obtain 

= -a) 

dl R' '' 


(3.2.19) 


where, q|r rc is in equilibrium with Cp at location R in the pellet. 

Therefore, assuming a linear isotherm for the gas adsorbed at the pore mouth, the 
following relationship between the adsorbate within the crystal and the gas phase in 
macro-pore is assumed: 

O.220) 


or. 




r R, 


KC. 


(3.2.21) 
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Here, K is the adsorption equilibrium constant/partitioning coefficient between macro- 
pores and crystal surface. 

The volume-averaged concentration {q) over the crystal is obtained by averaging 
equation (3,2. 1 7) and using equation (3.2.9) as follows: 

J=A+~BR^^ (3 2,22) 

5 

Similarly, the volume-averaged concentration over the pellet is obtained by averaging 
equation (3 .2. 1 6) and using Eq. (3.2.8); 

F,=C+-DR/ (32.23) 

5 

From equation (3.2. 1 6) the concentration at R ^Rp is evaluated as: 

C = (" + DR ^ (3,2.24) 


Combining equations (3,2.23) and (3.2.24) C and D are calculated as: 
5 — 

n = rc - C 1 
- 

and At _ / \ 

~ ^ 2 ^ 2 

Substituting the values of C and D in equation (3.2. 16) we obtain 


C = (— C 

^ V 2 P 




Differentiation of (3 .2.27) with respect to R we obtain 


(/; 

P 'I 


[C 

p 


The boundary condition at the pellet surface Eq, (3.2. 14) can be written as: 




= K„(C,-C 


P\R=R. 


Combining equations (3.2.28) and (3,2.29) we obtain 


p 

Rearranging equation (3.2.30) we obtain 


(3.2.25) 

(3.2.26) 

(3.2.27) 

(3.2.28) 

(3.2.29) 

(3.2.30) 
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(K .C , + ,) 

^ P 

(3.2.31) 

- 5 i) 


Substituting equations (3.2. 16) and (3.2.17) in equation (3.2.20) to obtain 


{A + = K{C + DR^) 

(3.2.32) 

Comparing term by term in equation (3.2.32), we obtain 


A = KC 

(3.2.33) 

BRl = KDR^ 

(3.2.34) 


Substituting the values of C and D from equations (3.2.26) and (3.2. 25) in equations 
(3.2.33) and (3.2.34), A and B are calculated as follows: 


A = K(-C„ p ) 

^9 P 9 ' 


(3.2.35) 


K 


B 


2R 


[C - C ]R 

^ P\R Ii„ ^ P J" 


(3.2.36) 


R^ 


On integration of equation (3.2. 1 1) we obtain an expression for q as follows: 




K [(’ + 



(3.2.37) 


Substituting the values of C and D from equations (3.2.25) and (3.2.26) we obtain an 
expression for q: 





3C 


p\R^R, 


] 


(3.2.38) 


Finally, substituting equation (3.2.38) and equation (3.2.21) in equation (3.2. 19), an 
expression for the mass balance of the average adsorbate phase is obtained as follows: 


3 q \ S 1) ^ ^ p. K ( ^ f! 

_= -yj-Ikc (8C, 

The equation is further rearranged by 
(3.2.3 1) in equation (3.2.39) to obtain 


- 3C 


p\R - R p 


) 


(3.2.39) 


substituting the value of Cp|R=Rp from equation 
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dq 


9D^KK, 




5D 

~T' 


(C.-CJ 


-) 


(3.2.40) 


The above equation is substituted in equation (3.2.15) to obtain a governing equation for 
the average pore concentration: 


dCp (C,-Cp)K„ 21DaK{l-a) 15D, 
R. 


Rc 


R: 


] = o 


(3.2.41) 


Similarly, on substitution equation (3.2.2) and equation (3.2.31) in equation (3.2.1), the 
equivalent simplified governing equation for the bed is obtained as follows: 


ac’t^r, eCt d. a^c, i5(i-E,)jc.£),(c,-c,) 


di 


8z 


Gj 




(3.2.42) 


The 2"'* order PDE equations (3.2.41) and (3.2.42) are the governing equations to predict 
the breakthrough characteristics of the bi-porous adsorbents like zeolite in a fixed packed 
bed. Thus, by assuming the volume-averaged quantities and parabolic concentration 
profiles in the macro and micro pores, the number of partial differential equation reduced 
3 to 2, which reduces the computational work/time to a large extent. 

3.3 Non-dimensionalisation of governing equations and Initial and Boundary 
conditions: 

The following dimensionless variables are used: 

C t 1 . . . 

Bulk gas concentration: C*fc Time: Distance in Axial 

^50 ^residence z 


Direction: = j » Radius of the pellet, R* 


R. 


Introducing the dimensionless variables in governing equation (3.2.42), the following non- 
dimensionalised equation is obtained: 


ac; 8c: , 1 a^c; 


dz Pe dz‘ 


n(c: -c:) 


(3.3.1) 
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where the constant terms are defined as 

1 ^ 15(1-6,)^ D/ 


r 


a 


5/:)„ 

^^RXi.K„+ ” 


R. 


) 


From the above expressions, it is easy to infer that Pe is nothing but the Peclet number, 
which compares convective flux to the flux by dispersion in the packed bed. Similarly r\ 
can be construed as the resistance to the overall mass transfer determined by film mass 
transfer coefficient and diffusivity in the macro-pore of the pellet. 

The non-dimensionalized boundary conditions are; 

at t*=0 Cj’ = 0 

at t>0 Z*=0 Cj = 1 


dc: 

dz* 

To non-dimensionlize equation (3.2.41), the following variables are used: 


— c. 

Average gas concentration inside the pellet, Cp = ; 

('bO 


Dimensionless time constant; 


t*= i 

The non-dimensionlized form obtained is as follows: 


dC 


f-+sic;-c:) = o 


(3.3.2) 


dt 

where the constant coefficient 5 is defined as 

15D, 


ffl p 


51) 


Rl 


r: 


A close inspection of the above coefficient 6 reveals that it is a product of two terms. 
While the first term includes the effect of film mass transfer coefficient and diffusivity 
within the pellet, the second term contains the adsorption equilibrium constant. It is 
important to point out that for a zeolite to be a sink for the adsorbing gas, the second term 
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must be negative. Thus, for a zeoilte-S02 system, the above term puts a constraint on, or 
in other words, defines a upper limit for the value of equilibrium constant K as: 

(3.3.3) 


K< 


\5l\Rl 


i?;27D,(l-a) 

For a given zeolite pellet and the adsorbing gas, the terms on the right hand side are 
known. Assuming Rp = 0.16 mm, Rc = 2x10"^ mm, De = 2.54x10"® mVs, Dc = 1x10"*^ 
mVs, and a = 0.64, K is evaluated to be 613. 

The non-dimensionlized initial and boundary conditions are as follows: 


0 


t*>0 R*=() 


= 0 


dC 


a/r 


0 


Equations (3.3.1) and (3.3.2) are the non-dimensionlized forms of the governing equations 
that are numerically solved in the present work to predict breakthrough times in a packed 
fixed bed consisting of bi-porous zeolite adsorbents. 

3.4 Numerical solution Technique: 

For the gas solid sorption reaction taking place we have two partial differential 
equations 3.3.1 and 3.3.2 containing two dependent variables C I and C* as a function of 

time and axial location. These equations are solved simultaneously by finite difference 
formulation using the NAG Fortran library. The Fortran codes are produced in the 
appendix. 
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Figure 3.1 ; Schematic diagram of sorption mechanism underlying non-porous model 


Partially reacted sub-grain 



Zeolite particle 



Figure 3.3 Schematic diagram of a bi-porous zeolite pellet 


CHAPTER 4 

EXPERIMENTAL STUDIES 


23 


Experiments were carried out under several operating conditions such as temperature, 
flow rate, particle diameter, amount of CaO and zeolite, and the inlet SO 2 concentration 
to study SO 2 sorption or adsorption. This chapter describes the experimental set-up and 
the experimental procedure. 

The experimental set-up used is the same as that designed and developed by 
Dasgupta (2001), except for a few modifications. A separate set-up was designed and 
developed to prepare precipitated CaO. 

4.1 Experimental Set up: 

The experimental set-up used in the process consisted of three sections as shown in 
Figure 4.1: a) gas mixing section, b) test section, and c) analytical section. In the gas 
mixing section two electronic mass flow controllers (Model PSFIC-I, Bronkhorst, 
Netherlands) were used to control and regulate the flow of SO 2 , O 2 or N 2 . A mixture of 
gases (SO 2 , O 2 or N 2 ) of the desired concentration was prepared by mixing the two 
individual streams in a gas-mixing chamber (135 mm long and 25 mm I.D.) made of 
stainless steel tube. The O 2 stream was purified using silica gel purifier to removal 
moisture in the gas. The gas-mixing chamber can be bypassed directly to the G.C for 
calibration as well as measurement of the inlet concentration of the gas mixture entering 
the reactor. 

The test section consists of a reactor and a furnace with a proportional 
temperature controller. An inconel reactor having an internal diameter 14mm and length 
150mm was used for the reaction of the sorbents with the inlet gas from the gas mixing 
section. The reactor was vertically mounted co-axially in a 200mm long tubular furnace 
The top end of the reactor was connected to the gas line using 1/4 swage lock fittings 
and specially mounted unions made of SS. Power to the fiimace was supplied through a 
15 A variac. The temperature of the furnace was controlled using a proportional 
temperature controller (Blue Bell, Kanpur). The temperature in the reactor was measured 
using a chromel-alumel thermocouple inserted in a S.S thermo well (O.D 3mm). The 
effluent gas stream from the reactor was passed to the analytical section. A manometer 





was also connected to the two ends of the reactor to measure the pressure drop across the 
sorbent bed. The analytical section consisted of a gas chromatography (GC) connected to 
a data station. A computer was connected to the data station to store the chromatogram 
and the peak areas. Thermal conductivity detector (TCD) was used to detect and measure 
SO 2 concentration in the gas stream. The calcination of the sorbent materials was done in 
a vertical quartz tube (50mm ID) having a purge arrangement. The calcination tube was 
vertically mounted in a movable fiirnace having diameter 35 cm and length 54cm. The 
temperature of the furnace was controlled by a proportional temperature controller 
(Model 961 HLTC, ionelectricals, Mumbai). The furnace temperature was measured 
using chromel-alumel thermocouple (Century Pvt. Ltd., Kanpur) of 12mm diameter and 
25 mm length. Power to the fiimace was supplied through a 15 A variac. 

Precipitation of CaCOn was done in a Pyrex reactor (35mm ID x 360 mm L). A PVC 
pipe with several small pore openings (-0.025 mm diameter) was used for bubbling CO 2 
in the aqueous slurry of Ca(OH) 2 . The flow of CO 2 was controlled by a needle valve and 
measured by a rotameter (0-5 slpm). Mixing of slurry was done with the help of a 
magnetic stirrer (Remi equipments, Kanpur) using Teflon needle. The precipitated 
solution was filtered using Whatman filter paper of 1 qualitative, 125mm-diameter 
circles. Filtered materials were dried in an oven (Mahendra Scientific Instrument Mfg., 
Kanpur). 

Activation of zeolites was done in a vacuum oven. Vacuum (650mm) was created 
with the help of a vacuum pump (Vacnik Rotary Vacuum pump. Model V2M5, Lawrence 
& Mayo (India Pvt. Ltd., Delhi). 

4.2 Experimental Procedure: 

All calcium based sorbents such as calcium carbonate (CaCOa), calcium hydroxide 
(Ca(OH) 2 ), calcium acetate(CaA), and calcium magnesium acetate (CMA) as well as 
molecular sieve zeolite (type 5A) were obtained from Merck Co, Germany. 

To obtain CaO form CaCOa, Ca(OH )2 and Ca(CHaCOOH )2 powders, calcination 
was carried out under inert N 2 environment with a flow rate of 1.6 slpm at 750, 850, 950, 
and 1050“C in a vertical, movable quartz furnace (refer Figure 4.2). The optimum 
calcination temperature for CaCOa and Ca(OH )2 was determined to be 950“G, while for 
the calcinations of Ca(CH 3 GOOH )2 a temperature of 1050°C was chosen. The 
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calcinations time in each case was 4 hours. The calcinations temperatures and time were 
determined after repeat of experiments (trial and error) so as to be consistent with 
stoichiomentric calculations for weight toss due to the evolution of CO 2 . 

The ppt-carbonate was prepared by carbonation-precipitation method in a slurry 
bubble column reactor (refer Figure 4.3). First, an aqueous slurry of calcium hydroxide 
was prepared by mixing a known quantity of CaO (2.5 gm) in approximately 130 cc of 
distilled water. The amount of water was used based on the reported solubility data for 
CaO in water [Guthrie et al.l901]. Calcium lignosulfonate was used as a surfactant 
medium in small quantity (2wt.% based on the weight of Ca(OH) 2 )) to uniformly disperse 
suspended calcium hydroxide in the slurry at 25“C held constant with the aid of a 
thermostat. Thereafter, carbonation was carried out by bubbling CO 2 at a steady flow rate 
of 1 00 cc/min for 4 hr through the Ca(OEr )2 slurry. The post-reaction slurry was filtered 
by #1 filter paper and dried in the oven at 75°C for 24 hr. In the case of CaCl 2 used as an 
additive in the sorbents, CaCl 2 was added to the filtered sample already dried for lOhr in 
the oven. Thereafter, the solid mixture was further dried for 14 hr before used for the 
sulfonation reaction. 

The zeolite materials obtained in the form of cylindrical pellets were crushed and 
sieved. The fractions investigated in this study comprised three size ranges; 1.115 mm, 
0.16mm, and 0.385 mm. N 2 was used as a carrier gas. Rest of the operating parameters 
were kept the same as used in the case of CaO 

The required quantity of the sorbents (either CaO or zeolite) was then fed to the 
reactor. The reactor was heated to the desired reaction temperature and further kept at 
that temperature for 2 hours so that the system was stabilized and a uniform temperature 
existed in the bed. The flow rate of the gas was adjusted using MFC. 

Mesh and quartz wools were put at both ends of the reactor tube so as to prevent 
any carry over of the sorbent. The thermo well for holding thermocouple was then placed 
co-axially in the center of the reactor tube. The gas line was connected using swage lock 
fittings. First, O 2 was allowed to pass through the reactor at a desired flow rate controlled 
by MFC. The identical flow rate reading by the rotameter installed at the downstream end 
of the tube, as that by MFC implied that the set up was free from any gas leakage. The 
sorbent was then heated to the required temperature. The concentration of the inlet 
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gaseous mixture in the line bypassing the reactor was measured by GC, prior to the start 
of the sorption reaction. The total experimental time included 2 hrs to pre-heat the bed to 
the desired temperature and 2 hours to stabilize. The reaction time varied between 2 and 6 
hrs depending upon the operating conditions including the type of sorbents, flow rates 
and the amount used. The concentration of exit gas from the reactor was measured by the 
GC. 

4.3 Product Analysis: 

The product gas was analyzed by gas chromatography using Porapak Q-column. The 
column (length 6-inch, I.D. 2.0mm) was used to separate O 2 and SO 2 using thermal 
conductivity detector (Model 5700.Nucon). Hydrogen at 22 ml/min was used as a carrier 
gas. The filament current was kept at 246 mA. The temperatures of the injector, detector 
and oven were maintained constant at 150°C, 160° C and 120° C, respectively. All the 
GC data were stored in a computer using a data station (Winchrom software) connected 
between the chromatograph and the computer. The areas obtained in the computer were 
in proportion to the concentrations as pre-defined in the software. After the run the 
product was analyzed for the amount of sulphate formed using UV spectrometer analysis, 
[.limited SEM analysis and the BET area measurements were also carried out for the 
samples of commercial calcium carbonate, ppt-calcium carbonate with and without 
CaCl 2 . 

4.4 Experimental variables: 

The experiments were performed in the packed bed reactor under different 
operating conditions such as bed temperature, inlet gas concentration, particle size, flow 
rate, and the amount and the type of the sorbents used. A comparative study was done by 
taking 2 g of the sorbent materials in powder form for each case. The gas concentration 
was kept constant at 10000 ppm, and flow rate at 600 slpm. To study the effect of 
reaction temperature on the removal of SO 2 the bed temperature was varied between 
750 ° C and 950 ° C, and the remaining parameters were kept constant. 

In case of zeolite, the reaction temperature was varied between 50°C and 100"C, 
The amount taken was 5 g. To study the effect of inlet gas concentration, particle size, 
and flow rate, one of the variables was changed and the remaining parameters were kept 
constant as shown in Table 4. 1 . 



Experimental 

variables 

Experimental conditions 

Gas flow rate 
(cc/min) 

300 

600 

Temperature 

("O 

50 

75 

100 

Sorbent particle 
radius 
(mm) 

0.164 

1.115 

0.385 

Inlet SO 2 
concentration 
(ppm) 

5000 

10000 

25000 


Table 4.1: DifTerent experimental variables for SO 2 sorption 





Figure 4.1: Schematic diagram of experimental set up 
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Figure 4.2: Arrangement for calcination at high temperature 
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Figure 4.3; Arrangement for Preparation of precipitated- CaO 
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CHAPTER 5 

RESULTS AND DISCUSSION 

In this chapter, both experimental and modeling results have been presented on the 
study of the sorption/adsorption behavior of SO2 on calcium oxide and zeolites sorbents In 
the case of CaO sorbents, analysis has been done for breakthrough as well as total sulphate 
conversion. 

5.1 Non-regenerative sorbent: 

The term breakthrough refers to the response of an initially clean bed (i.e. free of 
adsorbate) to an influent containing an adsorbate or a mixture of adsorbates. In the typical 
experiments performed, the reactor bed initially consisting of pure CaO is heated to a 
particular reaction temperature and fed with a gas containing mixture at a constant flow 
rate. During the initial stage of sorption reaction, most of SO2 entering the bed is removed 
by chemical reaction occurring at the surface of each sub-grain in the pores. During this 
stage, no or very low concentration of SO2 is observed at the reactor outlet. However, a 
product layer of CaS04 is gradually formed which offers resistance to the diffusion of SO2 
within the pores and preventing it from reacting with the CaO sub-grains. This resistance 
progressively increases till the diffusion of SO2 completely stops, implying that the overall 
sorption process is completed. Under this condition, the concentration of SO2 at the 
reactor outlet becomes equal to that at the inlet. In the present work, the transient change 
in SO2 concentration at the reactor outlet was determined using a gas chromatograph. 
Such transient behavior of concentration is called breakthrough and the corresponding 
concentration vs. time plot is called breakthrough curve. In the conversion analysis, the 
total amount of sulphate formed at the completion of the, calcium ion and any other rrietals 
present in the sorbent were determined by atomic absorption technique. 

For the model analysis the following quantities were calculated based on the 
experimental conditions; superficial gas velocity, residence time, bed porosity, density of 
the gas, Reynolds number, Schmidt number, dispersion co-efficient, Peclet number and 
mass transfer co-efficient. The corresponding calculations have been presented in the 
appendix. 





5.1.1 Effect of Temperature on sorbent: 

To determine the effect of reaction temperature on the SO 2 removal process, the 
experiments were earned out at 750, 850, and 950 °C. In all the runs, the flow rate was 
kept constant at 600 cc/min and 2 gm of the sorbents (powdered size) were used. The inlet 
SO 2 concentration was maintained constant at 10000 ppm. The breakthrough curves of 
CaO obtained from CaCOs, under various temperatures are described in Figure 5.1. As 
observed from the figure, the breakthrough times are approximately 25 min for both 
temperatures; 750 and 950 °C . However, the total or final reaction time (at the instance the 
reactor exit gas concentration asymptotically reaches the inlet gas concentration) is 
significantly greater at 950 (380 min) than at 750°C (150 min). The sulphate analysis of 
the reactant product was also carried out at the end of each run. As observed from Table 
5.1, the sulphate conversion increases from 24.63 % to 52.98 % as the temperature is 
increased from 750 to 950 °C. As also clearly observed from the tabulated results on 
breakthrough times and the total amount of sulfate formed, the sorption of SO 2 may be 
assumed to be insignificant at or below 750 “C. These results are consistent with those 
obtained by Dasgupta [2001]. In the present work, 950 °C was henceforth chosen as the 
reaction temperature for the sorption tests. As mentioned in Chapter 2, in an earlier study 
of Dasgupta [2001] sintering effects were observed to occur at high temperatures resulting 
in immediate breakthrough of the sorbents bed. In most of the experiments, the weight of 
the sample, gas flowrate and inlet concentrations were also henceforth kept the same (i.e. 
2 gm, 600 cc/min, and 10000 ppm, respectively) for a comparative study of various types 
of the sorbents. 

Figure 5.2 describes the breakthrough characteristics of CaO obtained from 
Ca(OH )2 at various reaction temperatures; 750, 850, and 950°C. The corresponding 
breakthrough times at 750, 850, and 950“C were approximately 10, 60, and 90 min, 
respectively. A comparative analysis suggests that CaO obtained from hydroxide 
performed significantly better than that obtained from carbonate under identical 
conditions. 

5.1.2 Sorbents obtained from acetate: 

CaO sorbents were also prepared by calcination of commercially available 
Ca(CH 3 COOH) 2 . The corresponding breakthrough characteristics of SO 2 at reaction 





temperatures ot 750, 850, and 950“C are described in Figure 5.3. As observed from the 
figure, the breakthrough time increased from 60 to 100 min with increase in the reaction 
temperature from 750 to 850®C. With further increase in the reaction temperature to 
950°C, the breakthrough time, however, decreased to 65 min. Referring to Table 5.1, the 
sulphate conversions at the reaction temperatures of 750, 850, and 950°C were determined 
to be 54.77%, 56.08% and 49.02%, respectively. The corresponding reaction times (for 
the reactor exit gas concentration equal to that at the inlet) were 235, 480, and 395 min, 
respectively. These results (breakthrough combined with sulphate analysis) clearly show 
that SO 2 sorption by calcium acetate based sorbents is most effective at 850°C, in 
comparison with 950 °C for carbonate and hydroxide based sorbents. 

Figure 5.4 describes the sorption efficiency of MgO obtained on calcinations 
from magnesium acetate at 750, 850, and 950°C. The breakthrough results clearly indicate 
that the sorption of SO 2 by MgO is insignificant, as the breakthrough occurs 
instantaneously, and in less than 10 min the reactor exit gas concentration reaches more 
than 60% of the inlet gas concentration. Furthermore, the breakthrough characteristics 
were observed to be invariant with the reaction temperatures. 

It is important to point out that Yiannis et al. (1993) carried out the sorption 
study on commercially obtained (from Cryogenic Technologies) calcium magnesium 
acetate (CMA), a mixture of calcium and magnesium acetates, and reported a significant 
superior performance of the sorbents in comparison with that of CaO obtained from the 
hydroxide or carbonate. In the present work, due to locally non-availability of this type of 
mixture or salt, a mixture of CaO and MgO obtained on calcination from their respective 
acetates was prepared in the laboratory by mixing in different proportions. Since the 
sorption capacity of pure MgO was observed to be insignificant, MgO was mostly used as 
an additive in the mixture of CaO and MgO. 

Figure 5.5 compares the breakthrough characteristics of various mixtures of MgO 
and CaO prepared in different ratios: 1:20, 1:8, 1:4, 1:2.6, and 8:1 those pure of MgO and 
CaO. A mixture containing higher amount of MgO than CaO was also tested for 
comparison purpose Table 5.2 summarizes the breakthrough times and total sulphate 
conversions corresponding to total reaction times for different mixture ratios. A 
comparison of the sorption characteristic of these CaO-MgO mixtures with that from CaO 
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obtained fiom calcium acetate (ref. Table 5.1) suggests that the mixing of MgO with CaO 
did not result in improvement of the breakthrough time or the sulphate conversion, though 
the overall performance of the mixture was observed to be marginally superior to that of 
carbonate based sorbents and to some extent, comparable with that of hydoxide based 
sorbents. This is in contrary to the observations or conclusion made by Yiannis et al. 
(1993) with regard to the superior performance of CMA in comparison with that od 
hydroxide or carbonate based calcium sorbents. 

A comparative performance of CaO sorbents obtained from different sources; namely: 
CaCOa, Ca(OH )2 and Ca(CH 3 COOH )2 is shown in Figure 5.6. The results are presented 
for the sorption carried out at the respective optimum reaction temperatures of these 
sorbents. As observed, acetate based sorbents exhibit longer breakthrough and reaction 
times than hydroxide based sorbents under identical operating condtions. On the other 
hand, the performance of carbonate based sorbents is inferior to hydroxide as well as 
acetate based sorbents. 

5.1.3 Precipitated (ppt) CaO based sorbents with and without CaCb (additives): 

The CaO-sorbents obtained from the calcination of the laboratory made ppt- 
CaCOs were tested for their sorption efficiency at reaction temperatures of 750 and 950 
°C. In some samples, additives CaCb were added in small quantities to determine the 
effect of moisture on SO 2 sorption. CaCl 2 being hygroscopic in nature, it was expected that 
some amount of moisture will be absorbed in the solid mixture. Figure 5.7 and 5.8 
compare the breakthrough characteristics of commercial and ppt-carbonate sorbents with 
and without CaCb mixed in various proportions, at 750 and 950 ®C, respectively. Tables 

5.3 and 5.4 compare the breakthrough and reaction times, along with the corresponding 
sulphate conversion for ppt-CaCOs with those for commercially available carbonate at the 
reaction temperatures of 750 and 950 °C, respectively. 

As observed from Table 5.3, the sulphate conversion at 750 °C for either 
sorbents (ppt and commercial) is barely 25%. However, addition of chloride increases the 
conversion to approximately 50% in 260 min. The corresponding breakthrough time 
increases from 35 min to 40 min. 

Similar to the sorption process at higher temperature, the total sulphate 
conversion increased significantly at 950 ^C. However, the performance of ppt-CaO was 



observed to be marginally higher than that of commercially obtained CaO, in terms of 
both reaction and breakthrough times. On the other hand, it was also observed that a 
certain amount of CaCla significantly improved the sorption capacity of ppt-CaO. As 
observed from Figure 5.8 and Table 5.4, the breakthrough time for the mixture of CaCl 2 
and CaO (0.05:1) considerably increased to as high as 120 min. The total reaction time 
was 410 min and the corresponding sulphate conversion obtained was 62%. Any other 
ratio higher or lower than 0.05 : 1 for CaCb and CaO yielded in shorter breakthrough time 
and lower conversion. 

Therefore, we conclude that in the present work, amongst various types of non- 
regenerative sorbents tested for SO 2 control under identical conditions, including 
commercial CaO obtained from carbonate, hydroxide, and acetate, commercial MgO, the 
solid mixture of CaO and MgO, laboratory prepared ppt-CaO with and without CaCb, ppt- 
CaO mixed with CaCb in 1 :0.05 ratio had the longest breakthrough and reaction times and 
highest sulphate conversion. Table 5.5 summarizes these comparative findings. 

5.1.4. Model results: 

As shown in the study of Dasgupta [2001], the model predictions assuming a 
porous structure of calcium based sorbents explained reasonably well the breakthrough 
characteristic of SO 2 . In the proposed mechanism, the overall SO 2 sorption by CaO in a 
fixed bed reactor under constant temperature conditions was determined by following rate- 
determining steps; 

a) diffusion from the bulk gas to the particle surface through a gas film siirrounding the 
particle (bulk diffusion). 

b) diffusion of the gas within the pores (pore diffusion) 

c) reaction of SO 2 at the surface of sub-grain or particle surface (chemical reaction) 

The diffusion of the gas from the bulk to the particle is dependent on the mass 
transfer co-efficient, which inturn is determined by the particle radius and gas flow rate. 
The chemical reaction rate constant has Arrehenius temperature depednecne. The 
diffusion of the gas in the pores is primarily dependent on the pore size of the sorbents. In 
line with the observations made by other investigators (Dogu[1982], Bruce[1989]) the 
pore size is assumed to vary during the sorption reaction, as the reaction products 
(sulphate) begin to gradually fill the pore, eventually resulting in pore-plugging. 
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Therefore, the eftective diffusivity in the pores also varies during the course of the 
reaction. The variation of the porosity as the reaction proceeds was assumed to be 

represented by an exponential function: a=af^e ^ where, ao= initial intraparticle 
porosity and d'and n are constants. Therefore, the proposed model had two adjustable 
parameters (<Jand n) to explain the experimental results on SO 2 breakthrough by the 
porous (mon-pore) sorbent particles. The rest of the parameters such as effective pore 
diffusivity taking into account Knudsen diffusivity, gas molecular diffusivity, gas film 
mass transfer coefficient and reaction rate constants were calculated based on the initial 
existing conditions. 

The model predictions of SO 2 breakthrough on various types of the sorbents 
under varying operating conditions are shown in Figures 5.7 and 5.8. As observed, the 
predictions compare reasonably well with the experimental data. The breakthrough 
characteristics of various calcium based sorbents is further explained by analyzing model 
predicted variation in the porosity, diffusivity, and gas concentration within the pores of 
the CaO particles during the course of the sulfonation reaction. In the simulation, the 
molecular diffusivity of SO 2 was calculated be 6.39 XIO’^ and 8.36X10'^ at 950°C and 
750“C, respectively. The reaction rate constant was calculated to be 5 and 8.46 m/s at 750 
and 950'’C, respectively and the initial porosity of the CaO particle was assumed to be 
0.52 [Dogu, 1980]. Based on the intitial porosity and pore size, initial Knudesn (4.58x10 '* 
m^/s) and effective diffusivity (3.33x10'^ mVs) in CaO sorbents were also calculated. The 
results are shown in Appendix). 

As the sulfation proceeds the porosity of the CaO particles gradually decreases at 
the reaction temperatures. Incomplete conversion occurs when the porosity within the 
grain asymptotically becomes zero. Figure 5.9 describes the variation in porosity at 750°C, 
as predicted by the model. As also observed, the porosity decreases at a slower rate for p- 
CaCOs mixed with CaCb in a mixing ratio 0. 1/1.0, compared with the sorbents of other 
mixing ratios. The decrease in the porosity is highest for c-CaO. Figure 5.10 describes the 
model predicted corresponding decrease in the effective diffusivity of the reacting gas in 
the pores of CaO particle at temperature 750”C. The porosity and hence difbisivity were 
gradual for 0.1/1 than other mixing ratios. Figure 5.1 1 describes the average concentration 



profile in the pore of CaO particles, as predicted by the model. As observed from the 
figure, concentration profiles decrease slowly for mixing ratio, 0.1:1 that supports the 
higher sulphate conversion. Similar observations were observed at 950°C for the case of 
mixing ratio 0.05:1.0. The average concentration profiles in the pore are shown in Figure 
5.13. Figures 5.14 and 15 depict the variations in the porosity and diffixsivity, respectively. 
The behavior observed for the mixture ratio 0.05:1 at 950°C are identical to those for the 
mixture ratio of 0.1:1 .0 at 750°C. 

Since in this work, the BET area of the sorbent was measured to be low (2~10 
m“/g) the combined experimental and model results suggest that sorption mechanism of 
SO 2 on CaO based sorbents is mainly pore diffusion controlled. 

5.2 Regenerative sorbents (zeolite) 

5.2.1 Adsorption temperatures for zeolites: 

Similar to the experiments carried out on calcium based sorbents to determine optimum 
temperature for maximum SO 2 sorption, the tests were also carried out on 5A zeolites at 
various temperatures to obtain the optimum temperature for maximum adsorption of SO 2 . 
Figure 5.15 describes the breakthrough characteristics of SO 2 during dynamic adsorption 
by zeolite materials in a fixed packed bed at reaction temperatures of 100, 75, and 50“C. 
The gas flow rate was kept constant at 300 cc/min and the gas inlet SO 2 concentration was 
maintained at 10000 ppm. The zeolites materials commercially available in the form of 
pellets were crushed to prepare a 5 gm test sample having average particle size 1.115 mm. 
It was observed that at temperatures greater than 100°C and lower than 50°C, there was 
practically no significant adsorption, resulting in instantaneous (< 3 min) breakthrough of 
the bed. At temperature as high as 200®C, zeolite pellets were, in fact, burnt off. 

As observed from Figure 5.15, the optimum temperature for SO 2 removal was 
obtained at 75‘’C. The breakthrough time at this temperature was approximately 50 min, 
while it took more than 90 min for the reactor exit gas concentration to approximately 
reach the inlet concentration (10000 ppm). On the other hand, at 100 °C, the breakthrough 
time decreased to 35 min. The corresponding adsorption time was 80 min. At 50 °C, the 
breakthrough of the bed occurred in less than 10 min. The adsorption temperature for all 
the experiments hence after was selected as 75 "C. The suitability of 5 A zeolites within a 
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temperature range of 50 to 100®C for the maximum adsorption of SO 2 has also been 
reported el sewhere T sibranska [1999]. 

Model predictions were done under identical experimental conditions. As 
described in the section on theoretical analysis, the mathematical model is based on two 
governing equations; one for the bulk gas phase and the other in the macro-pore. The 
concentration profiles within the zeolite crystal due to the radial diffusion in the intra (or 
micro)-pores of crystals was averaged and the partitioning between the interface of the 
crystal surface and marco-pores was assumed to be determined by a parameter, called 
equilibrium constant. This dimensionless equilibrium constant has been used as an 
adjustable parameter in the model, wherever necessary, to explain the experimental data. 
The remaining variables such as gas film mass transfer coefficient around the pellet, 
dispersion coefficient in the packed bed and macro-pore diffusivity incorporating Knudsen 
effects and intra-particle porosity were calculated. The temperature effects on variation in 
gas, pore or crystal diffusivities were incorporated, as necessary. The macro-pore size for 
5A zeolites was assumed to be 0.3 pm based on reported value Do et al [1998]. The 
sample numerical calculations are reported in Appendix. 

Under the experimental conditions chosen for the results on Figure 5.15, the 
gas film mass transfer coefficient and the effective pore diffusivity were calculated to be 
0.05 m/s and 2.54x10"^ mVs, respectively. As observed, a reasonable good agreement is 
observed between the model predictions and the data within the experimental and 
computational errors. The values for partition coefficient K at 50, 75 and 100 were 
adjusted at 45, 13, and 4, respectively. The decrease in the values of K with increase in the 
temperature is also consistent with an independent work Kopac et al [1997]. It may be 
pointed out that in the aforementioned study, the equilibrium constant was defined 
between the bulk gas phase and crystal surface concentrations, as against between the 
macro-pore and the crystal surface concentrations in this work. Assuming Arrhenius type 
temperature dependence, the activation energy for adsorption/partitioning was calculated 
to be 46.4 KJ/mol. The model predicted values of the average concentration within the 
macro-pores of the pellets at these temperatures are described in Figure 5.16. The 
characteristics of these curves are consistent with the proposed mechanism. The macro- 



pores get saturated with the adsorbate at the instance the tube exit concentration reaches 
the inlet gas concentration. 

5.1.2 Effect of particle size: 

The experiments were carried out for three different particle (pellet) sizes: 0.160 
mm, 0.385 mm and 1.115 mm under varying gas flow rates and inlet SO 2 concentrations. 
The test samples having a particle size of 0.16 mm were in nearly powdered form, while 
those having higher size particles were in spherical pellet forms. The amount of zeolite 
materials was taken as 5 g in each case. In one test run, the inlet SO 2 gas concentration 
and flow rate were kept constant at 5000 ppm and 600 cc/min, respectively, and in the 
second run, these values were 10000 ppm and 300 cc/min, respectively. Figure 5.17 and 
5.18 describe the effects of particle sizes on the breakthrough characteristics for these 
runs. As observed from the figures, decrease in particle size results in increase in 
breakthrough and also, total adsorption times. However, as the particle size is decreased 
(below 0.385 mm), the reductions in these times are observed to be minimal. From these 
experimental data it was concluded that lower the particle size, superior is the performance 
of zeolites, in terms of breakthrough time or adsorption capacity, the remaining operating 
conditions being identical. Thus, if pressure drop or agglomeration of the particles is not a 
limiting factor, the adsorbents in powdered form should be used in the packed bed. 

Model predictions were done to explain the effect of particle size on the 
observed breakthrough times. It was found that as per the constraint imposed by eq. 
(3.3.3) on the upper limit of K for the term 5 to be negative, the equilibrium constant value 
(model parameter) was adjusted to approximately 613 corresponding to the particle size of 
0.16 mm. The gas film mass transfer co-efficient was calculated to be 0.05 m/s for 1.115 
mm size (radius) particle, as compared to 0.23 m/s for 0. 16 mm size particle. 

The superior performance of relatively smaller size zeolite materials can be 
explained as follows. With the increase in mass transfer co-efficient due to smaller particle 
radius, the amount of gas or the flux of SO 2 (mol/m ^ -s) diffusing into the adsorbents 
increases. This, in turn, results in higher adsorption rate. As the particle radius is 
decreased for the same amount of sorbent the interfacial area per unit volume (a*) also 
increases. Hence, there is a combined effect of mass transfer rate as well as the interfacial 
area per unit volume on the total SO 2 flux entering the pores of the pellet. The increase in 



mass transfer coefficient and the interfacial area due to reduction in particle size, however, 
could not explain the corresponding improved breakthrough characteristics. 

Referring to the assumptions in the model analysis, fhe crystal radius of zeolite 
materials was assumed to be constant (2xl0'^m). Since there was no temperature variation 
between two conditions, the gas diffusivities (both bulk as well as Knudsen) were 
unchanged. Since intra-particle porosity (determined by bulk and solid densities) may be 
assumed to be constant for both particle sizes, the number of crystals is less in smaller size 
particle than in larger size particle, for the same quantity of zeolites. Therefore, for the 
same amount of the gas in the macro-pore as well as within crystal, the concentration 
(moles/cc) within the crystal volume, q, is more. The partioning coefficient/equilibrium 
constant between the gas concentration in the macro-pores and within the crystals was 
accordingly increased from 14 to 605 to explain the experimental observations, as the 
pellet size decreased from 1.115 to 0.16 mm. The model predicted values of the average 
concentration within the macro-pores of the pellets at these pellets size are described in 
Figure 5.20. 

5.2.3 Effect of flow rate: 

To observe the effect of flow rate the experiments were carried out at optimum 
selected temperature i.e. 75°C on 5 g of adsorbents having particle size 0.16mm. The inlet 
gas concentration was kept constant at 5000 ppm. Figure 5.21 describes the breakthrough 
characteristics at these operating conditions. As observed from Figure 5.21, the 
breakthrough time decreases from 98 to 28 min as the gas flowrate increases from 300 to 
600 cc/min. This decrease is attributed to an early saturation of the bed due to relatively 
larger amount of SO 2 at higher flowrate. In this case, the data could be explained by the 
model with a marginal change in the model parameter K (from 605 to 601 i.e. less than 
1 %). 

5.2.4 Effect of concentration; 

To observe the effect of concentration the experiments were carried out for 
different gas inlet concentrations: 10000 and 5000 ppm, the remaining operating variables 
were kept the same as for the conditions of the previous section. Figure 5 .21 describes the 
breakthrough characteristics at these concentrations. As observed from /fjie figure, the 



breakthrough time increases from 12 min to 28 min with decrease in the inlet gas 
concentration from 10000 ppm to 5000 ppm. 

This can be explained as follows. With decrease in the inlet gas concentration for 
identical flowrates, the total amount of SO 2 entering the maco-pores of the pellet is less. 
Therefore, saturation of the pores occurs in relatively longer time. In this case also, the 
data could be explained by the model with a marginal change in the model parameter K 
(from 605 to 598 i.e. less than 1 %). The model predicted values of the average 
concentration within the macro-pores of the pellets at these concentrations are described in 
Figure 5.23. 

5.2.5 Effect of adsorbent amounts: 

To observe the effect of weight of the adsorbing materials (5 A zeolite), the 
experiments were carried out on two different amounts; 2g and 5 g. The particle size in 
each case was 0.16 mm. The inlet gas concentration and flow rate were kept constant at 
5000 ppm and 600 cc/min, respectively. 

Figure 5.24 describes the SO 2 breakthrough charcteristics for these amounts. 
As observed, the breakthrough time increases from 8 min to 28 min, with increase in the 
weight of zeolites from 2 to 5 g. The increase in the breakthrough time is obviously due to 
approximately 45% increase in the total surface area for adsorption of SO 2 . Once again, 
the data could be explained by the model with a marginal change the parameter K (from 
605 to 610 i.e. less than 1%). 

Table 5.8 summarizes the various input parameters used in the model, such as 
difflisivity, mass transfer coefficient, etc corresponding to the various experimental 
operating conditions used in this work, such as gas flowrate, temperature and 
concentrations, etc. The last column lists the adjusted values of the model parameter K for 
these conditions to explain the observed breakthrough characteristics of SO 2 on 5A 
zeolites. 

5.3 Pressure drop across the packed bed: 

Referring to Figure 4.1 describing the schematic of the experimental setup, a U- 
type manometer (CCI 4 as a medium) was used to measure the pressure drop across the 
adsorbent bed for various sorbent amounts, particle sizes, and the gas flow rates. The 
measurements are reported in Tables 5.9 and 5.10 for particle size 1.115 mm and 0.16 


mm, respectively. As observed from Table 5.9, the measured pressure drop across the 
packed bed increases with increase in the flow rate. For the particle size 1.115 mm the 
pressure drop across the bed at flow rate 300 cc/min was apprpximately 7.04 N/m^, which 
increases to 83.63 N/m at the flow rate 600 cc/min. The effect of particles size on the 
pressure-diop can also be observed from the tables. The pressure drop across the bed 
increases from 7.04 N/m to 83.63 N/m ^ with decrease in the particle size from 1.115 mm 
to 0.16 mm. 

The tables also compare the experimentally measured pressure drop with those 
theoretically obtained using well known Ergun equation for the pressure drop in a packed 
fixed bed [McCabe et al, 1997], The detailed calculations are described in Appendix. A 
significant difference between the experimental and calculated values may be due to large 
pressure drop across the mesh (wool) installed at the reactor inlet to support the 
adsorbents, and at the inlet to prevent any carryover of the powdered particles with the 
flowing gas. 

h'or industrial purposes, adsorbents of larger pellet sizes may be desirable to 
minimize the pressure drop across the bed. The pressure drop calculation is also required 
to determine the gas blower capacity. 

5.4 Comparative performance of calcium based sorbents and zeolites: 

A comparative performance of calcium based sorbents and 5A zeolites for control 
of SO 2 is described in Figures 5.25 and 5.26. Comparison has been made between these 
materials (non-regenerative and regenerative types) at their respective optimum reaction 
temperatures (i.e. for zeolites, T = 75 ®C and for CaO, T= 850 °C). However, the weight of 
the sorbents/adsorbents and the gas flowrates were kept the same for comparison. Figure 
5.25 compares the breakthrough characteristics of 5 A zeolites with those of CaO obtained 
from acetate for various gas flowrates. It may be recalled that acetate based sorbents were 
found to be superior in comparison with those obtained from carbonate or hydroxide. 
Figure 5.26 compares the breakthrough characteristics of these materials for varying inlet 
gas concentrations, As observed from each figure, calcium acetate may be considered far 
superior to zeolites, as the breakthrough and total reaction times for the former are at least 
10 or in some cases, 100 times greater than those for the latter. However, due to their 
inherent regenerative characteristics, zeolites may still be preferred to non-regenerative 



types of sorbents such as CaO. Use of a two-bed system, similar to PSA (pressure swing 
adsorption), may also be advantageous in the case of zeolite adsorbents. 

5.5 Model simulation of control of a large-scale SO 2 emission: 

Model simulations were carried out to determine the total amount of 
adsorbents in a fixed bed to control a large-scale SO 2 emission from an industry. With the 
help of the model, the SO 2 breakthrough and reactions times were also predicted. The 
simulations were done for various gas flow rates between 3000 and 10000 slpm. A vertical 
reactor (adsorbent bed) of 0.3 m I.D. and 3.15 m height containing 120 kg of zeolites 
(average particle size 1.115 mm) was assumed for the control of SO 2 emission. The 
emission level was assumed to be 10000 ppm. The bed was operated isothermally at 
75° C. Figure 5.27 describes the simulated breakthrough curves for varying flow rates of 
the gaseous effluents. As observed from the figure, the breakthrough times were predicted 
to be approximately 220, 90 and 10 min corresponding to the gas flowrates of 3000, 5000 
and 10000 slpm, respectively. Thus, with the help of the above-simulated results, 
changeover or regeneration time to put a bed of fresh adsorbent can be estimated. For 
example, if the gas flow rate is 3000 slpm, a changeover time of approximately 4 hours 
may be required. Hence, the simulated results are useful in designing an industrially sized 
adsorber. 
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Figure 5.1: Effect of temperature on SO2 sorption by CaO obtained from carbonate 
(Qgo = 6 cc/min, Qq = 600 cc/min, Cg, = 10000 PP™’ W = 2 g) 
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Figure 5.2: Effect of temperature on SO 2 sorption by CaO obtained from hydroxide 
(QgQ = 6 cc/min, Qq = 600 cc/min, Cg, = 10000 ppm, W = 2 g) 
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Figure 5.4: Effect of temperature on SO2 sorption by MgO obtained from acetate 
(Qso = 6 cc/min, Qq = 600 cc/min, Cg, = 10000 ppm, W = 2 g) 
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Figure 5.13: Model predictions of variation of diffusivity with time for additive (CaCl2) in ppt-CaO 

(Qso = 6 cc/min, Qq = 600 cc/min, T = 950®C, W = 2 g) 
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Figure 5.14: Model prediction for average concentration profiles in the pore on SO2 sorption 

(Qgo = 6cc/min, Qo»= cc/min, T = 950®C, W = 2 g ) 
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Figure 5.15; Effect of temperature on SO2 adsorption by zeolite 
(Rp = 1.115 mm, Q = 300 cc/min, Cgjj^j . = 10000 ppm, W = 5 g ) 
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Figure 5.17: Effect of pellet size on SO2 adsorption by zeolite 
( Q = 300 cc/min, W = 5 g, T = 75®C, Cgjjjjg^ = 5000 ppm, De = 2.54e-6 m^/s, Dc = le-12 m 
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Figure 5.18: Effect of pellet size on SO2 adsorption by zeolite 
( Q = 300 cc/min, W = 5 g, T = 75®C, Cgi^iet = 1000® = 2.54e-6 m^/s, Dc = le-12 m^s) 
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Figure26: Comprasion of SO 2 sorption by zeolite and CaO 
(zeolite: T = 75®C, CaO: T = 850®C, W = 5 g) 





S.No. 

Sample 

Temperature 

("Q 

Break- 

through 

Time 

(min) 

Reaction 

Time 

(min) 

Sulphate 

conversion 

(%) 

1 

CaCOj 

750 

25 

150 

24.63 

2 

CaCOs 

950 

25 

380 

52.98 

3 

Ca(OH )2 

950 

90 

400 

48.60 

4 

Calcium 

acetate 

750 

60 

235 

54.77 

5 

Calcium 

acetate 

850 

100 

480 


6 

Calcium 

acetate 

950 

65 

395 

49.02 


TableS.l: Breakthrough time, reaction time and sulphate conversion 
of different commercial sorbents at different temperatures 


S.No. 

CaO + MgO 
(g) 

Temperature 

rc) 

Break- 
throug 
h time 
(min) 

Reaction 

Time 

(min) 

Sulphate 

Conversion 

(%) 

1 

0.25 + 1.75 

850 

20 

140 

18.55 

B 

1.25 + 0.75 

850 

60 

410 

70.23 

B 

1.5 + 0.5 

850 

40 

240 

40.13 

B 

1.75 + 0.25 

850 

65 

1 

370 

64.69 


Table5.2: Breakthrough time, reaction time and sulphate conversion 
of different mixer (CaO + MgO) ratios of sorbents 








































S.No. 

Types 

of 

sorbents 

Break 

through 

Time 

(min) 

Reaction 

time 

(min) 

Final 

wt. 

(g) 

Sulphate 

conversion 

(%) 

1 

c-CaCOs 

25 

150 

2.618 

24.63 


Ppt-CaCOs 

35 

165 

3.42 

28.79 

3 

Ppt-CaC 03 + 
CaCl2 
(0.05/1.0) 

■ 

40 

220 

3.314 

41.07 

4 


40 

260 

3.735 

54.68 

5 


40 

65 

2.847 

26.65 


Table 5.3: Breakthrough time, reaction time and sulphate conversion 
of different sorbents at temperature TSO^C 















S.No. 

Types 

of 

sorbents 

Break 

through 

Time 

(min) 

Reaction 

time 

(min) 

Final 

wt. 

(g) 

Sulphate 

conversion 

(%) 

1 

c-CaCOs 

25 

380 

4.332 

52.97 


ppt-CaCOs 

30 

440 

4.557 

56.19 

3 

Ppt-CaC 03 + 
CaCl 2 
(O.Ol/l.O) 

55 

290 

3.812 

45.69 

4 

ppt-CaC 03 + 
CaCb 
(0.02/1.0) 

105 

400 

3.9.13 

44.61 

5 

ppt-CaC 03 + 
CaCl2 
(0.05/1.0) 

120 

410 

3.992 

62.18 

6 

ppt-CaC 03 + 
CaCl 2 
(0. 1/1.0) 

50 

370 

1 

3.914 

50.49 

7 

ppt-CaCOs + 
CaCl 2 
(0.15/1.0) 

50 

360 

3.991 

38.01 


Table 5.4: Breakthrough time, reaction time and sulphate conversion 
of different sorbents at temperature 950‘*C 














S.No. 

Types 

of 

sorbents 

Temperature 

rc) 

Break- 

through 

Time 

(min) 

Reaction 

Time 

(min) 

Sulphate 

Conversion 

(%) 

1 

Calcium 

carbonate 

950 

25 

380 

52.98 

2 

Calcium 

hydroxide 

950 

90 

400 

48.6 

3 

Calcium 

acetate 

850 

100 

480 

56.08 

4 

CaO + MgO 
(1.25+0.75) 

850 

60 

410 

70.23 

■ 

ppt-CaCOa 

950 

30 

440 

56.19 

6 


950 

120 

410 

62.18 


Table 5.5: Comparative performance of various non-regenerative sorbents 




















S.No. 

Temperature 

CC) 

Break- 

through 

Time 

(min) 

Adsorption 

Time 

(min) 

1 

50 

10 

40 

2 

75 

50 

90 

3 

100 

35 

80 


Table 5.6: Breakthrough and adsorption time at different 
temperature 


S.No. 

Particle 

size 

(mm) 

Cone. 

(ppm) 

Flow rate 
(cc/min) 

Break- 
through 
time (min) 

Adsorption 

time 

(min) 

1 

1.115 

10000 

300 

50 

90 

2 

1.115 

5000 

300 

92 

170 

3 

0.385 

10000 

300 

55 

105 

4 

0.385 

5000 

300 

100 

170 

5 

0.160 

10000 

300 

50 

125 

6 

0.160 

5000 



300 

[ 

95 

180 


Table 5.7 Breakthrough time and adsorption time at different flow 
rate and concentration for various particle size 

















































S.No 

Pellet 

radius 

(mm) 

Wt 

(g) 

Temp 

("C) 

Flow 

rate 

(cc/min) 

Cone, 

(ppm) 

Molecular 

diifusivity 

(Dm) 

(m*/s) 

Effective 

diffusivity 

(De) 

(m*/s) 

Ciystal 

diffusivity 

(Dc) 

(mVs) 

Mass 

transfer 

coefficient 

(Km) 

(m/s) 

Partitioning 

coefficient 

K 

■ 

1.115 



600 

5000 

2.07e-5 

2.54e-6 

le-12 

0.07 

13 

■ 



50 



1.85e-5 

2.3e-6 

2.57e-13 

0.05 

45 

■ 

1.115 


75 

300 

10000 

2.07e-5 

2.54e-6 

le-12 

0.05 

13 

■ 

1.115 

1 

100 

n 



2.77e-6 

3.24e-12 

0.05 

4 

■ 

0.16 

a 

75 

600 

5000 

2.07e-5 

2.54e-6 

le-12 

0.23 

610 

■ 


5 

75 

600 

5000 

2.07e-5 

2.54e-6 

le-12 

0.23 

605 

7 



75 

300 

5000 

2.07e-5 

2.54e-6 

le-12 

0.20 

601 

8 

0J6 

5 

75 

600 

10000 

2.07e-5 

2.54e-6 

le-12 

0.23 

598 


Table 5.8: Various Input parameters for model predictions 
(Rc = 2x10'’ m, P = 1 atm, a = 0.64) 






















































S.No. 

Wt. 

(g) 

■ 

Flow 

rate 

(cx:/min.) 

Mano- 

meter 

Height 

(CCI 4 ) 

(cm) 

Vbeci“V2»o 

Vw 

(Sb) 

Data 

pressure 

P=pgh 

(N/m^) 

Pressure 
(Ergun eq) 
at (Sb) 
(N/m^) 

Error 

(Pexp-P*e)*10() 

P.xp 

(%) 

1 

2 

3 

300 

4.0 

0.354 

624.1 

7.04 

98.87 

D 

2 

3 

600 

10.0 

0.354 

1560.2 

14.09 

99.09 

B 

5 

Hi 

300 

8.7 

0.192 

1357.2 

138.22 

89.82 

B 

5 

B 

600 

19.5 

0.192 

3042.0 

276.44 

90.91 

B 

B 

B 

300 

9.5 

0.152 

1482.5 

409.13 

72.40 

6 

7 

B 

600 

20.0 

0.152 

3120.1 

818.25 

73.77 


Table 5.9: The experimentally and theoretically calculated pressure drop 
across the reactor for 1.115mm particle 


S.No. 

Wt. 

(g) 

Height 

(h) 

(cm) 

Flow rate 
(cc/min) 

Mano- 

meter 

Height 

(CCI 4 ) 

(cm) 

Data for 
Pressure 
P=pgh 

'^bed"V2eo 

Vw 

(^b) 

Pressure 
(Ergun eq) 
(NW) 

Error 

(P=xp-Pthe)*100 

Pexp 

(%) 

1 

B 

4 

300 

27.0 

4212.4 

0.5149 

83.63 

98.01 

H 

B 

■ 4 

400 

65.0 

10141.0 

0.5149 

111.51 

98.90 

B 

B 

B 

450 

88.0 

13729.4 

0.5149 

112.45 

99.17 

B 

B 

8 

150 

75.0 

" 1 

11701.2 

0.3938 

291.93 

97.50 

5 

B 

8 

200 

83.5 

13105.3 

0.3938 

389.25 

97.03 


Table 5.10: The experimentally and theoretically calculated pressure drop 
across the reactor for 0.16mm particle 

































































CHAPTER 6 
CONCLUSIONS 


In the present work, a combined experimental and theoretical study was done to 
understand the kinetics and mechanism of SO 2 removal by 5A zeolite materials. The 
emphasis was on determining the suitability of the commercially available zeolites in the 
adsorption of SO 2 under varying operating conditions. Experimental screening of various 
non-regenerative types of Ca-based sorbents was also carried out to determine their 
comparative sorption performances measured in the breakthrough time and total sulphate 
conversion. Most of the sorbents, such as CaCOa, Ca(OH) 2 , Ca(CH 3 COOH )2 were 
commercially obtained. Some of the sorbents were prepared or modified in the 
laboratory. 

The main findings of this work can be summarized as follows. CaO obtained on 
calcination of calcium acetate was found to be superior amongst all sorbents used in this 
work, with relatively lower reaction temperature (850°C), longer breakthrough time (100 
min), longer reaction times (480 min) and also, higher sulphate conversion (~70%). CaO 
obtained on calcination of CO 2 precipitated carbonate in the laboratory exhibited 
marginally superior performance relative to that obtained commercially. However, the 
addition of small quantity of CaCl 2 in ppt-CaCOa (0.05 .1) significantly improved the 
performance, as the breakthrough time increased from 25 min to 120 min and the 
corresponding sulphate conversion from 50% to 60%. CaCb being hygroscopic in nature, 
the incorporation of trace amount of moisture in the sorbents resulted in decrease of the 
volume of the sulphate layers. The pore closure phenomenon was thus delayed. 
Magnesium based sorbents were found to be the least suitable for SO 2 sorption. 

In the case of adsorption by zeolites, the various experimental conditions included 
the reaction temperature, inlet gas concentration, flow rate, particle size and the amount 
of zeolites. The optimum adsorption temperature was determined to be 75®C. The 
adsorbents of smaller size exhibited relatively longer breakthrough and reaction times. A 
mathematical model for zeolites (bi-porous) was developed to predict the SO 2 
breakthrough characteristics in a packed bed under a wide range of operating conditions. 
The theoretical analysis incorporated the effects of gas diffusion in the adsorbent bed, in 



the macro-pore and within the crystal. A simplified approach based on linear driving 
approximation and average radial concentrations significantly reduced the computation 
time. The model had only one adjustable parameter, the equilibrium or partitioning 
constant or coefficient. The model predictions were found to be in good agreement with 
the experimental data. 

6.1 Scope of future study; 

Further search on or development of novel sorbents is required. As the sorption 
process suffers from premature termination due to pore closure mechanism, an optimum 
gas solid contacting pattern is required, as in a fluidized bed. Additionally, simulated 
moving bed, high gravity contactor can also be tried to improve upon the mass transfer 
rate from bulk of the gas to the surface of the particle. 

In this work, the adsorption of SO 2 by 5 A type zeolite was carried out. A study on 
novel adsorbents, such as different types of zeolites such asl3X molecular sieve, 
molecular sieve carbon (MSC), and activated carbon fibers (ACF) is required. The 
regenerative characteristic of the adsorbents should be investigated so as to reuse the 
adsorbent and make the process continuous. A surface analysis of the sorbent materials 
before and after the reaction, using SEM, XPS should be carried out to determine the 
change in the morphology/structure of the particles during the course of the reaction. 
Similarly, pore size distribution and BET area measurements should also be done to 
determine change in the pore structure of the particles. 
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APPENDIX 


The following calculations are used in the model based on the experimental conditions; 
gas flow rate = 600 cc/min 


superficial gas velocity = {Flcwrate){I{^Temp) 

i22A){6Q)(Pressure){nd\U) 


(Al) 


where, 

flow rate is in lit/min 
R = universal gas constant in J/kg K 
Temp = temperature in K 
Pressure m N fm^ 

Cross sectional area in 


superficial velocity = 


(0.6)(8.314)(348)(4) 
(22.4)(60)(10 1 325)(;r(0.0 14)' ) 


0.0828 m/s 


Residence time of gas in the reactor = 


length 

vel 


Residence time = ^ = 0.9661 s 

0.0828 

Concentration ~ ^ 

(R){Temperature) 

^ 001322 ( 00264 )^^,, 224 

(8.314)(348) 

Bed porosity 

Amount of sorbent = 5 gm 
density of sorbent = 670 kg/ m' 

S X 1 0~^ 

Volume of CaO = — — = 7.4626x10-® m' 

670 


TV 


Volume of the packed bed = -(0.014)' (0.08) = 12.308 x 10'® m' 

4 


„ , . (12.308- 7.4626) xlO-® . 

Bed porosity e. = = 0.3936 

^ ^ 12.308x10-® 


(A2) 


(A3) 


(A4) 



Bed porosity = 


(12.308- 7.4626) xlO~^ 
12.308x10“' 


0.3936 


Intra-particle porosity 


(— — L)/ 

d„ d 



( 


1 


1.97 


)/ 


I 

0.71 


(A4) 

(A5) 


= 0.64 
Density of gas 

_ (^pressnre)(mol.wt) 
(Temperature)(R) 


(101325)(64xl0~^) 

(8.314)(348) 


= 2.2413 kg/m^ 


Reynolds number based on particle diameter 

(density)(sup vel)(particlediameter) 

JtvC — 

(viscosity) 


(A6) 


(A7) 


2.2413x0.0828x1.6x10-^x2 

2x10“' 


2.96 


Schmidt number 


Sc = 


(viscosity) 

(density)(molecidardiffusivity) 


2x10“' 

2.2413 x 2.073x10“' 


0.43 


Dispersion co-efficient 

The dispersion co-efficient (Dz)\s expressed as 


dispersion - coefficient 


(par tide. dia)(vis cos ity) 
(density)(Pe) 


(A8) 


(A9) 


where Pe is Peclet Number. 

The Peclet number are correlated with Reynold’s number and Schmidt number by the 
correlation 


1 _ 0.3 0.5 

^"Rex.S'c l+3.8/(RexJ?c) 


(AlO) 



0.3 


Pe 


(2.96)(0.43) 


05 

(1 + 3.8/(2.96x0.43)) 


Pe = 2.77 


Z)z= (2xl.ll5xl0“^)x2xl0-' ..,. 92 , 

^ = 7.2xlO®mV5 

2.2413x2.77 

For 0.008 < Re < 400 

0.28 < Sc <2.2 

The mass transfer co-efficient is calculated using the following correlations; 
6’/j = 2.0 + l.l(6’c)°''(Re)°^ 


_ (2 + 1 . !(&)'’'' (Re)° " ){mol.diffusmty) 

— 

{part. diameter) 

^ _ (2 + 1. 1(0.43)°^^ (2.96)°^)(2.073xl0~^) 
” ~ (2x1.6x10-^) 

K„, = 0.233 m/s 


(All) 


Calculation of diffusivity: 

The effective diffusivity inside the pores is given by a combination of both Knudsen and 
molecular diffusion. 

The Knudsen diffusivity is given by 




2xr 

3 ^tvM’ 




(A12) 


where 

Rpore== pore radius in m 
M = Molecular weight 
T = Temperature in K 
Dk = Knudsen diffusion 

The Molecular diffiisivity for a binary gas mixture is given by the Chapman - Enskog 
equation given by 


10'^[1.084-0.249(— + 


D = 

■ m 


M. M, 


1 i ^ 
)2]7’2( 


1 


1 




fCT 


(A13) 



Here 


Dm = Molecular diffusivity mVs 
r = T emperature in K 

Ma, Mb = Molecular weight of the two interacting species 
P = Total Pressure NW 

KT 

/( )= Collision integral, a fiinction oi kT! s where k is Boltzaman constant. 

e 

•^ = 335 . 4 ^: 

K 


^ = 71 AK 
K 


^ = V335.4x71.4 
K 

= 2.248 K 
KT 

/(— )= 0.46 
€ 


-^_A±Ib 


‘AB 


’SO^-N^ 


0.4112 + 0.3798 


nm 


= 0.3955 nm 

Substitute all the values in A(12) 


10-^[1.084-0.249( 


1 1 i,J,l 1 


D. 


64 28 


) 2 ] 72 (. 


64 28 


(101325)(0.3955)"(0.46) 

=2.073x10-^ mVs 

The combined diffusivity contributed by both molecular and Knudsen diffiision is given 
by 
1 


_ = _L _L 

D ~ D 

A m 


(A14) 


D = Diffusion co-efficient 



The actual diffusion path in the pores will not be equal to the distance in the radial 
direction, but will be quite tortuous. The actual path depends on the pore structure and the 
tortuosity factor is defined as the ratio between the actual diffusion path length and the 
net distance in the direction of flux or radial distance. 

The effective difflisivity inside the pores is found by using the formula 



where 

a = intraparticle porosity 
D = difflisivity co-efficient 
r = tortuosity factor 
Calculation of diffusivities 
For macro-pore radius 0.3 pm 

The Knudsen difflisivity = (97)(0. 3x10'^) 

= 6.78x10'^ mVs 

The Molecular difflisivity = 2.073x10'^ m^/s 

The Diffusion co-efFicient is given by 

1 ^ 1 1 
D ~ 6.78x10-' 2.07x10-' 



D = 1.58x10 '/w^/sec 

The effective difflisivity is calculated as follows 



r 


where r = 4, a = 0.64 
1.58x10-' xO.64 


= 2.54x10^/m"/s 


Calculation for activation energy; 


K = K.^e 


RT 


A(16) 



In -^ = for Ki = 45 at 323 K 

and K2= 13 at 348 K 
E (activation energy ) was calculated 

,45 E , \ 1 , 

13 8.314 323 348 


E = 46.4 KJ/mol 

The temperature dependence of Crystal dififusivity is Arrhenius-type expression 

A(17) 

at temperature 348 K the crystal difflisivity is assumed to be 1x10'*^ mVs 
Hd = 12.09 Kcal/mol Ref: [C.Wu et al., 1984] 

12.09x4,18x10^ 1 1^ 

(A) 32 ,r =2-57xlO-”».Vi 


Sample calculation for the conversion analysis: 

After the reaction is over, the sample was analysed for the amount of sulphate present in 
%. The sorbent after reaction has un-reacted CaO as well as CaS04 formed after the 

reaction. Suppose the analysis revealed a 30.8 % <S'04 ” after the reaction. 

100 gm of sample gives 53.41 gm S04'' 


3.992x53.41 

3.992 gm (final wt.) will give = loo 


2.132 


Now 96 gm of sulphate gives 136 gm of CaS04 

2.132 gm will give gm of CaS04 = _ 3 0205 

96 

Reaction CaO ~ ^ CaS04 

56 gm CaO gives 136gm CaS04 

2gm ► 136x2 = 4.857 

56 

4.857 gm CaS04 will obtain at 100% conversion 



3.0205 gm CaS 04 will obtain at 


100x3,0205 

4.857 


62.188 


Pressure Drop Calculations: 

Experimental Pressure drop across the reactor 150 mm long and 14mm ID 

AP = pgish 

where 

p = density of manometric liquid 
g = acceleration due to gravity 
Ish - difference in the level of manometric liquid 
The reading in the manometer is 4 mm of CCU 
AP= (4.0e-2)(9.81)(l. 592x10^) 

= 624.06 N/m' 

Pressure drop across sorbent bed as per Ergun equation (McCabe et al, 1960) 
^_ 150//i;(l-e,)^I 

where 

p = viscosity of the gas 
Vz = superficial velocity 
Ej, = bed porosity 

L = length of the packed bed 
ij>^ = sphericity (assumed = 1) 

Dp = diameter of particle 

AP = 0 50)(2e - 5)(0.0414)(1 - 0.354)^(3e - 2) ^ ^ 

(2xl.l5e-3)'(0.354)' 

Reynold’s Number based on tube diameter 

Re = 0.4xl0-)(0 . 0 1 H) g,211 3 ) ^ ^ 

(2c -5) 

The pressure drop across the tube 


(A18) 


(A19) 



(A20) 


Df 

where 

// = viscosity 
Vz = superficial velocity 
1 = length of the reactor 
Dt = diameter of the tube 

^ _ (32)(2e-5)(0.0414)(146g-: 
(I4e-3f 

erimental ^^bed ^^reactor 


= 0.0 197 AT //w^ 



C Last change: I 13 Feb 2002 7:25 pm 
C Last change: SA 27 Feb 2001 8:48 pm 

IMPLICIT DOUBLE PRECISION (A-H.O-Z) 
IMPLICIT INTEGER (l-N) 


C Parameter statements for principal dimensions 

PARAMETER ( NPAR = 10 , NSPAR = 1 , NSTEQ = 2) 
PARAMETER ( LRW = 5000000 , LIW = 70000 ) 

PARAMETER ( NPTMAX = 100001 ) 

C Dimension statements and Externals 

DIMENSION Y(NSTEQ, NPTMAX), YPRIME(NSTEQ.NPTMAX) 
DIMENSION INFO(18),RWORK(LRW).IWORK(LIW) 

DIMENSION RPAR(NPAR),IPAR(NSPAR) 

DIMENSION ZGRID(NPTMAX) 

EXTERNAL DIFF,FUNC,BNRY,ESUB 

COMMON/prop/poros,dl,vel,dlen,rp,dka,pden,bdkm,dporos1 ,pdiff1 , 
1 rc,cdiff,beta 


OPEN(UNIT=1 .FILE=’*c:\agupta1\code\i.txt",STATUS=’'UNKNOWN’’) 
OPEN(UNIT=2.FILE="c:\agupta1 \code\2a.txt", STATUS="UNKNOWN'’) 
OPEN(UNIT=3,FILE="c:\agupta1\code\3a.txt",STATUS=”UNKNOWN") 
OPEN(UNiT=4,FILE="c:\agupta1\code\4a.txt",STATUS=”UNKNOWN") 

READ(1,*) Gflow.gvis,tdia,dlen,Pmolin,gmw, Temp, Pres, bdm.rp, 

1 pden,dka,wt,pdiff1,cdiff,rc,Tactual,npts,tgrad,dporos1 

wRITE(4,*) Gflow,gvis,tdia ,dlen . Pmolin,gmw,T emp,Pres,bdm,rp, 

1 pden,dka,wt,pdiff1 ,cdiff,rc,Tactual,npts,tgrad.dporos1 
CLOSE(I) 

PI=4.0D0*ATAN(1.0D0) 

Rcons=8.314 [universal Gas constant (J/gmol-K) 

Temp=Temp+273.0 IK 

Vel=4.0*Gflow*Rcons*Temp/(22.4*60.0*Pres*PI*tdia**2.0) Isuperficial gas Velocity, m/s 

Gden=Pres*Gmw*1 .0e-3/(rcons*Temp) !kg/m3 

Re=gden*vel*rp*2.0/gvis 

sc=Gvis/(gden*bdm) 

pe=1/(0.3/(re*sc)+0.5/(1+3.8/(re*sc)))!peclet no from Yang 
DI=(2.0*rp*gvis)/(Gden*pe) ! m2/s 
bdkm=(2.0+(1 .1 *(sc)**0.33*(re)**0.6))*bdm/(2.0*rp) 

Cgin=pres*Pmolin*gmw*1 .0e-3/(rcons*temp) I kg/m3 
WRITE(4,*)"Gden,re,sc,pe,DI,bdkm” 

WRITE(4,*)Gden,re,sc,pe,DI,bdkm 
Time=dlen/Vel !sec 

WRITE(4,*)"time,vel.cgin" 

WRITE(4,*)time,vel,cgin 



poros=(dlen*pi*tdia**2/4.0-wt*1e-3/pden)*4.0/(cllen*pi*tdia**2) 

WRITE(4.*)"poros=".poros 

c1= (1 .0-dporos1)*27*cdiff*dka/(rc**2) 

d1= 15*pdiff1/rp**2 

beta=bdkm*rp*’*2/((bdkm+(5*pdiff1/rp))*pdiff1*dporos1)*(c1-d1) 

WRITE(4,*)" beta=-’. beta 

WRITE{4,’*)"c1=",c1 

WRlTE(4,*)"d1=".d1 

IF(beta.Ge.0)WRITE(6,*)"beta is POSITIVE" 

RTOL =1.00-5 
ATOL =1.00-5 

C NPTS = number of points for axial grid (= 1 + intervals) 

C RTOL = relative tolerance for POESAC; an absolute tolerance is 

C calculated later 

C Preliminaries are done; now set things up for POESAC 
ICORO = 0.0 

C Integration times and initial values for states 
C at the discretization points ZGRIO(I) 


T=Q 

ntout=1 00000 


C Initial values for state variables 
C Impurity concn; use an exponential function in z 
C Pressure: use a linear profile for p(z) 


OO 30 1=1, NPTS 

ZGRIO(I) = OBLE(l-1)/OBLE(NPTS-1) 
Y(1. 0=0.0 
y(2, 0=0.0 

30 CONTINUE 


ZGRIO(NPTS)= 1.000 
Y(1,NPTS) = 0.0 
y(2.npts) = 0.0 


C Prepare for nagOOl ; Insert the INFO0 values. 

C ' Indicate that we have a banded Jacobian: info(6). 

C Ask POESAC to evaluate the initial yprime: info(1 1). 

C All others are "default" settings 


OO 40 J = 1,18 
INFO{J) = 0 



40 CONTINUE 
INFO(6) =1 


C Supply the bandwidths of the iteration matrix. 

IWORK(1) = 2*NSTEQ - 1 
IWORK(2) = 2*NSTEQ - 1 

C Integration loop. Call the implicit integrator PDESAC. 

WRITE(2,1 02)tOUT*time/60.0,y(1 ,1).y(1 .npts/2),y(1 .npts) 
WRITE(3.102)tounime/60.0,y(2.1),y(2,npts/2).y(2.npts) 

tout=tgrad 

DO 110 IOUT= 1.NTOUT 
c WRITE(2.*)iout,"1" 

IF(tout*time.ge.tactual)GO to 1 1 

Tout=1.05*tout 

50 CALL nag001(T.ZGRlD.Y,YPRIME,TOUT.INFO,RTOL,ATOL. 

1 IDID,RWORK,LRW,lWORK,LIW.RPAR,IPAR, 

2 NSTEQ.NPTS.ICORD.DIFF.FUNC, 

3 BNRY.ESUB) 


C Write results regardless of the integration status, 
c Stop if IDID is negative, 
c write(2,*)iout 

IF (IDID .Lt. 0) THEN 
WRITE(1. 160) IDID 
GO TO 140 
ENDIF 

WRITE(2,1 02)tOUT*time/60.0,y(1 .1).y(1 .npts/2).y(1 ,npts) 
WRITE(3,102)tounime/60.0,y(2,1),y(2.npts/2).y(2,npts) 

110 CONTINUE 

1 1 continue 

101 FORMAT(f15.5.(F15.9)) 

102 FORMAT(f15.5.3(F15.9)) 

140 WRITE(1 ,150)IWORK(1 1),IWORK(12).IWORK(13).IWORK(14) 
150 FORMAT(/1x;Number of steps taken so far :=M5/ 

1 IX/Numberof function calls :-,l5/ 

2 lX.'Numberof Jacobian calls :=',I5/ 

3 IX/Number of error test fails :-,l5/) 

160 FORMAT(/1X, 'Integration Failed with IDID =’,15) 


C End of main program 



c STOP 
END 


SUBROUTINE FUNC(T,X,U.UX.DUXX.FVAL.NSTEQ.RPAR.1PAR,IRES) 
IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

IMPLICIT INTEGER (l-N) 

DIMENSION U(1),UX(1).DUXX(NSTEQ.1).FVAL(1),RPAR(1).IPAR(1) 
COMMON/prop/poros,cll,vel,dlen,rp,dlka,pder>,bdkm,dporos 1 ,pdiff 1 , 

1 rc,cdiff,beta 

alpha=-1 -O/poros 
theta=dl/(vel*dlen*poros) 

eta =1 5*(1 .0-poros)*bdkm*pdiff1*dlen/(poros*rp**2*vel* 

+ (bdkm+(5*pdiff1/rp))) 

Fval(1)=alpha*ux(1)+theta*duxx(1.1)-eta*(U(1)-U(2)) 


Fval(2)=-beta*(u(1)-u(2)) 

RETURN 

END 


SUBROUTINE DIFF(T,X,U.DVAL,NSTEQ,RPAR,IPAR,IRES) 
IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

IMPLICIT INTEGER (l-N) 

DIMENSION U(1).DVAL(NSTEQ,1).RPAR(1),IPAR(1) 


C Supply the multiplier of second derivative. 

DVAL(1,1) = 1.0 
Dval(2,1) = 0.0 

RETURN 

END 


SUBROUTINE BNRY(T,X.U.UX, DELTA, NSTEQ.RPAR.IPAR.IRES.IBNRY) 
IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

IMPLICIT INTEGER (l-N) 

DIMENSION U(1).UX(1),DELTA(1),RPAR(1),IPAR(1) 

IF (IBNRY .EQ. 1) THEN III Left boundary (Z=0) 

DELTA(1)=U(1)- 1.0 
delta(2)=u(2)-1.0 

ENDIF 


IF (IBNRY .EQ. 2) THEN 


!!! Right boundary (Z=1) 



DELTA(1) = UX(1)- 

0.0 

delta(2) = ux(2)- 

0.0 

ENDIF 

RETURN 

END 

SUBROUTINE 

ESUB 

RETURN 

END 


Input file for the bi~porous model: 

0.3 2.0e-5 0.014 0.08 0.005 64.0 75 101325 2.073e-5 0.16e-3 
670 601 5.0 2.54e-6 le-12 2e-6 12000 100 0.01 0.64 

Gflow, gvis, tdia, dlen, pmolin, gmw. Temp, Pres, bdm, rp, 
pden, dka, wt, pdiffl, cdifF, re, Tactual, npts, tgrad, dporosl 


Determination of various parameters in so2vipp.txt 

1. flow in std litre per minute 

2. viscosity(kg/m-sec) 

3. tube dia in m 

4. tube lenghth in metre 

5. mol fraction 

6. molecular weight (g/gmol) 

7. temperature in "C 

8. pressure in Wm^ 

O.molecular diffusivity (m^/sec) 

10. radius of partical (m) 

1 1 .bulk density of zeolite (kg/m^) 

12. adsorption equilbirium constant 

13. porosity of bed 

14. diffusion of SO 2 in pores (mVsec) 

15. diffusion of SO 2 in crystal (m^/sec) 

16. crystal radius (m) 

17. actual time (sec) 

18. number of axial grids 

19. time gradint 

20. intraparticle porosity 



